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PREFACE 



This publication describes the IBM 
Systein/360 Time Sharing System IBM FORTRAN 
IV language referred to in this manual as 
FORTRAN IV. A reader should have some 
knowledge of an existing FORTRAN language 
before using this publication. The publi- 
cation FORTRAN General Information ^ Form 
F28-8074,, is a useful source for such 
knowledge. 

The material in the FORTRAN IV publica- 
tion is arranged to provide a quick defini- 
tion and syntactical reference to the var- 
ious elements of the language by means of a 
box format. In addition, sufficient text 
describing each element,, with appropriate 
examples of possible use, is given. 

There are four appendixes which give 
additional information useful in writing a 
FORTRAN IV source program. They are: 

A: Table of Source Program Characters 
B: Other FORTRAN Statements Accepted 

by IBM FORTRAN IV 
C: FORTRAN Supplied Subprograms 
D: Sample Programs 



First Edition 

Significant changes or additions to the specifications contained in this 
publication will be reported in subsequent revisions or Technical 
Newsletters . 

This publication was prepared for production using an IBM computer to 
update the text and to control the page and line format. Page 
impressions for photo-offset printing were obtained from an IBM 1403 
Printer using a special print chain. 

Requests for copies of IBM publications should be made to your IBM 
representative or to the IBM branch office servina vour locality. 

A form is provided at the back of this publication for reader's 
comments. If the form has been removed, comments may be addressed to 
the IBM Corporation, Time Sharing System/360 Programming Publications, 
Department 504, 2651 Strang Blvd., Yorktown Heights, N. Y. 10598 

© International Business Machines Corporation 1966 



CONTENTS 



INTRODUCTION 5 

The IBM System/360 Time Sharing System 

FORTRAN IV- 5 

Features of the Time Sharing System 

FORTRAN IV. ..-....», 5 

ELEMENTS OF THE LANGUAGE ........ 7 

Statements ........ - . ..... 7 

Coding FORTRAN Statements - Card 

Input. 7 

Coding FORTRAN Statements - 

Keyboard Input ........... 8 

Initial Lines- .......... 9 

Continuation Lines ........ 9 

Constants. ...-10 

Integer Constants .......... 10 

Real Constants. .11 

Complex Constant, -.-,-..... 12 

Logical Constants .13 

Literal Constants - - .13 

Variables. , 13 

Variable Names- .,...,,.... lU 
Variable Types and Length 

Specifications ........ . , . lU 

Type Declaration by the Prcidefined 

Specification - . 15 

Type Declaration by the IMPLICIT 

Specification Statement. 15 

Type Declaration by Explicit 

Specification Statements ...... 16 

Arrays ... ...... 16 

Declaring the Size of an Array. ... 18 

Arrangement of Arrays in Storage. . . 19 

Expressions. 19 

Arithmetic Expressions- ....... 19 

Arithmetic Operators . 20 

Logical Expressions ......... . 23 

Relational Operators 24 

Logical Operators- ... 24 

ARITHMETIC AND LOGICAL ASSIGNMENT 

STATEMENT ...-..»-....... 27 

CONTROL STATEMENTS .... . ., - . . . . 29 

The GO TO Statements - . . , ,., ..... 29 

Unconditional GO TO Statement. . . 29 

Computed GO TO Statement ..... 30 

The ASSIGN and Assigned GO TO 

Statements- ..,...,..-.. 31 

Additional Control Statements,, - - . - - 32 

Arithmetic IF Statement. ..... 32 

Logical IF Statement .«.-.-- 33 

DO Statement .... - 34 



CONTINUE Statement 37 

PAUSE Statement. ......... 38 

STOP Statement .38 

END Statement 39 

INPUT/OUTPUT STATEMENTS. ......... UO 

Basic Input/Output Statements. ..... 40 

READ Statement 41 

The Form READ (a,x). 42 

The Form READ {a,b) List ..... 44 

The Form READ (a) List 4 5 

Indexing I/O Lists ......... 46 

Reading Format Statements. .... 47 

WRITE Statement .47 

The Form WRITE (a,x) 48 

The Form WRITE (a,b) List 49 

The Form WRITE (a) List.. ..... 50 

FORMAT Statement. . . . . 50 

G Format Code. 52 

Numeric Format Codes (Ii,F,,E,D) . . 56 

I Format Code. . . . . .57 

F Format Code 58 

D and E Format Codes 58 

L Format Code 59 

A Format Code 59 

Literal Data in a Format 

Statement - 61 

H Format Code. ........... 62 

X Format Code 63 

T Format Code. . . . . .64 

Scale Factor - P . . . 64 

Carriage Control 66 

Additional Input/Output Statements ... 67 

END FILE Statement .67 

REWIND Statement 67 

BACKSPACE Statement. 67 

SPECIFICATION STATEMENTS 68 

The Type Statements 68 

IMPLICIT Statement 68 

Explicit Specification 

Statements. 70 

Adjustable Dimensions. 72 

Additional Specification Statements. . , 73 

DIMENSION Statement 73 

COMMON Statement 74 

Blank and Labeled Common ..... 76 

EQUIVALENCE Statement 77 

SUBPROGRAMS. 80 

Naming Subprograms ,80 

Functions .80 

Function Definition. ....... 81 

Function Reference 81 

Statement Functions .81 

FUNCTION Subprograms 83 

Type Specification of the 
FUNCTION Subprogram 84 



RETURN and END Statements in a 

Function Subprogram 85 

SUBROUTINE Subprograms .86 

CALL Statement , 87 

RETURN Statement in a SUBROUTINE 

Subprogram- ... ^ ....... 88 

Multiple ENTRY into a Subprogram . 89 

Additional Rules for Using ENTRY - 91 

The EXTERNAL Statement 92 

FORTRAN Supplied Subprograms, .... 93 

BLOCK DATA Subprogram 93 

APPENDIX A: SOURCE PROGRAM CHARACTERS . 95 

APPENDIX B: OTHER FORTRAN FEATURES 

ACCEPTED BY FORTRAN IV. 96 

READ Statement .......... 96 

PUNCH Statement- 96 

PRINT Statement. . 97 

DATA Initialization Statement. . . 97 

DOUBLE PRECISION Statement .... 98 



Arguments of a FUNCTION or 
SUBROUTINE Program Enclosed by 
Slashes 98 

APPENDIX C: FORTRAN SUPPLIED 
SUBPROGRAMS 99 

Mathematical Function Subprograms. ... 99 

Machine Indicator Tests .102 

The EXIT, DUMP, and PDUMP Subprograms. .102 

EXIT Subprogram 102 

DUMP Subprogram. 103 

PDUMP Subprogram 103 

APPENDIX D: SAMPLE PROGRAMS ,104 

Sample Program 1 .104 

Sample Program 2 .105 

INDEX ,. .112 



ILLUSTRATIONS 



FIGURES 



TABLES 



Figure 1. FORTRAN Coding Form 8 

Figure 2. Sample Program 1 104 

Figure 3. Sample Program 2 107 



Table 1. Insurance Premium Codes. ... 18 
Table 2. Determining the Mode of an 

Expression Containing Variables of 

Different Types and Lengths ...... 21 

Table 3. Valid Combinations With 

Respect to the Arithmetic 22 

Table 4. Mathematical Function. 

Subprograms ..... 99 



INTRODUCTION 



THE IBM SYSTEM/ '360 TIME SHARING SYSTEM FORTRAN IV 

The IBM Time Sharing System/360 FORTRAN IV is comprised of a 
language, a library of subprograms* and a compiler. 

The FORTRAN language is especially useful in writing programs both in 
conversational and nonconversational mode for scientific and engineering 
applications that involve mathematical computations. In fact,, the name 
of the language - FORTRAN - is derived from its primary use; FORmula 
TRANslating. 

Source programs written in the FORTRAN language consist of a set of 
statements constructed from the elements of the language described in 
this publication- 

The FORTRAN compiler analyzes the source program statements and 
transforms them into an object program that is suitable for execution on 
the IBM System/360. In addition,, when the FORTRAN compiler detects 
errors in the source program, appropriate error messages are produced. 
At the user's option a complete listing of the source program is 
produced. 

The FORTRAN compiler operates under control of Time Sharing 
System/360, which provides the FORTRAN compiler with input/output and 
other services. Object programs generated by the FORTRAN compiler also 
operate under System/360 Operating System control and depend on it for 
similar services. 

The IBM Time Sharing System/360 FORTRAN IV language is compatible 
with and encompasses the American Standards Association (ASA) FORTRAN, 
including its mathematical subroutine provisions. 



FEATURES OF THE TIME SHARING SYSTEM FORTRAN IV 

The Time Sharing System/360 FORTRAN IV is a further development of 
previously implemented FORTRAN systems and contains many of the features 
of these systems. In addition, the following features facilitate the 
writing of source programs and reduce the possibility of coding errors: 

1. Variable Attribute Control: The attributes of variables and arrays 
may now be explicitly specified in the source program. This 
facility is provided by a single explicit specification statement 
which allows a programmer to: 

a. Specify storage length. 

b. Explicitly type a variable as integer, real, complex,, or 
logical* 

c. Specify the dimension of arrays. 

d. Specify data initialization values for variables. 
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2, Adjustable Array Dimensions: The dimensions of an array in a 
subprogram may be specified as variables; when the subprogram is 
called, the absolute array dimensions are substituted. 

3. Additional Format Code: An additional format code - G - can be used 
to specify the format of numeric and logical data. Previously 
implemented format codes are also permitted. 

U, Mixed Mode: Expressions may consist of constants and variables, of 
the same and/or different types and lengths, 

5. Named I/O List: Formatting of input/output data is facilitated by 
reading and writing operations, without reference to a FORMAT 
statement or list. 

6. Spacing Format Code: The T format code allows input/output data to 
be transferred, beginning at any specified position. 

7. Literal Format Code: Apostrophes may be used to enclose literal 
data. 



ELEMENTS OF THE LANGUAGE 



STATEMENTS 



Source programs consist of a set of statements from which the 
compiler generates machine instructions,, constants, and storage areas- 
A given FORTRAN statement effectively performs one of three functions : 



Causes certain operations to be 
branch) . 



performed (e.g, , add,, multiply. 



2- Specifies the nature of the data being handled. 

3i Specifies the characteristics of the source program. 

FORTRAN statements are usually composed of certain F0RT31AN key words 

used in conjunction with the basic elements of the language: constants, 

variables,, and expressions. The five categories of FORTRAN statements 
are as follows: 

1. Arithmetic and Logical Assignment statements : Upon execution of an 
arithmetic or logical assignment statement,, the result of calcula- 
tions performed or conditions tested replaces the current value of 
a designated variable or subscripted variable. 

2. Control State ments ; These statements enable the user to govern the 
flow and terminate the execution of the object program. 

3. Input/Output Statements ; These statements,, in addition to control- 
ling input/output (I/O) devices, enable the user to transfer data 
between internal storage and an I/O medium. 

4. Specification Statements : These statements are used to declare the 
properties of variables, arrays,, and subprograms (such as type and 
amount of storage reserved) and to describe the format of data on 
input or output. 

5. Subprogram Statements ; These statements enable the user to name and 
define functions and subroutines. 

The basic elements of the language are discussed in this section. 
The actual FORTRAN statements in which these elements are used are 
discussed in following sections. 



CODING FORTRAN STATEMENTS - CARD INPUT 



The Statements of a FORTRAN source program can be written on a 
standard FORTRAN coding form,. Form X28-7327 (Figure 1) . FORTRAN 
statements are written one to a line from columns 7 through 72,. If a 
statement is too long for one line* it may be continued on as many as 19 
successive lines by placing any character, other than a blank or zero,, 
in column six of each continuation line. For the first line of a 
statement;, column six must be blank or zero* 

Columns 1 through 5 of the first line of a statement may contain a 
statement number consisting of from 1 through 5 decimal digits. Leading 
zeros in a statement number are ignored. The statement numbers may be 



Elements of the Language 



assigned in any order; the value of statement numbers does not affect 
the order in which the statements are executed in a FORTRAN program. 

Columns 73 through 80 are not significant to the FORTRAN compiler and 
may, therefore,, be used for program identification,, sequencing, or any 
other purpose. 

Comments to explain the program may be written in columns 2 through 
80 of a line, if the letter C is placed in column one. Comments may 
appear anywhere within the source program. They are not processed by 
the FORTRAN compiler, but are printed on the source program listing. 

Blanks may be inserted where desired to improve readability. 
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Figure 1.. FORTRAN Coding Form 



CODING FORTRAN STATEMENTS - KEYBOARD INPUT 



It is desirable to free a conversational keyboard operator from 
strict positional requirements when typing in a FORTRAN source program. 
The following conventions for statement numbers, text starting posi- 



tions, and continuation lines are accepted when input is from an on-line 
keyboard. 



Initial Lines 

If a line is the initial line of a statement, it may have a statement 
number. (The statement number,, if any, must appear on the first line of 
the statement- ) The numeric statement number must be the first nonblank 
material in the line. It can start in any column, and is terminated 
after five adjacent columns, or by the occurrence of a nonblank, 
nonnumeric character,, whichever happens first. 

If a statement has a statement number,, the text of the statement 
begins with the first nonblank character following the statement number, 
unless this character is a horizontal tab. If a tab is used to separate 
the statement number from the text, the text begins with the first 
nonblank character following the tab. 

If a statement does not have a statement number, the text of the 
statement begins with the first nonblank, nonnumeric character of the 
line* unless this character is a horizontal tab. 

If a tab is used to begin the line, the text starts with the first 
nonblank character following the tab. 



Continuation Lines 

A line of input is a continuation line, rather than the initial line 
of a statement, if the last character (blanks included) of the last 
preceding noncoroment line was a '-* (EBCDIC 60). 

A continuation line in keyboard input may not have a statement 
number. The text of the line begins with the first character (blank or 
not) of the line, unless this character is a horizontal tab. If a tab 
is used to begin the line, the text starts with the first character 
(blank or not) following the tab. 

Caution is needed in the use of '-' at the end of a line, and tab (in 
alphameric constants) or the letter C, at the beginning of a line, to 
avoid conflict between the FORTRAN text and the continuation, comment, 
and tab conventions of keyboard input. 



Elements of the Language 



CONSTANTS 

A constant is a fixed, unvarying quantity. There are three classes 
of constants — those that deal with numbers (numerical constants), 
those that deal with truth values (logical constants),, and those that 
deal with literal data (literal constants). 

Numerical constants may be integer, real, or complex numbers; logical 
constants may be .TRUE. or .FALSE.; literal constants may be a string 
of alphameric and/or special characters enclosed by quotes. 



INTEGER CONSTANTS 



r ■ — 1 

I Definition | 

^ ^ 

Integer Constant - a whole number written without a decimal point. 
It occupies four locations of storage- 
Maximum Magnitude: 2147U83647, i.e., (23^-1). 

L J 



An integer constant may be positive, zero, or negative; if unsigned, 
it is assumed to be positive. Its magnitude must not be greater than 
the maximum and may not contain embedded commas. 

Examples ; 

Valid Integer Constants: 



91 

173 

-2147483647 

-12 

Invalid Integer Constants : 

0.0 (contains a decimal point) 

27. (contains a decimal point) 

3145903612 (exceeds the allowable range) 

5,396 (embedded comma) 



10 



REAL CONSTANTS 



r 1 

I Definition | 

|. _ . ^ 

Real Constant ; - a number with a decimal point optionally followed 
by a decimal exponent, or an integer constant followed by a decimal 
exponent. This exponent may be written as the letter E or D 
followed by a signed or unsigned, one- or two-digit integer 
constant. A real constant may assume one of two forms: 

1. From one through seven decimal digits optionally followed by an 
E decimal exponent. This form occupies 4 storage loccitions. 

2. Either one through seven decimal digits followed by a D decimal 
exponent or 8 to 16 decimal digits optionally followed by a D 
decimal exponent. This form occupies eight storage locations 
and is sometimes referred to as a double precision constant. 

Magnitude: (either form) or IG'^^a through 16^^ (i.e., approxi- 
mately 10^5) . 



A real constant may be positive, zero, or negative (if unsigned, it 
is assumed to be positive) and must be of the allowable magnitude. It 
may not contain embedded commas. The decimal exponent E or D permits 
the expression of a real constant as the product of a real constant 
times 10 raised to a desired power. 

Examples ; 

Valid Real Constants (4 storage locations): 

+ 0- 

-999.9999 

0.0 

5764.1 

7-OE+O (i.e., 7.0 x 10° = 7.0) 

19761. 25E+1 (i.e., 19761-25 x 10^ = 197612.5) 

7E3 

7,E3 

7,0E3 > (i.e., 7.0 x 10^ = 7000.0) 

7.0E03 

7.0E+03 

7.0E-03 (i.e., 7.0 x 10-^ = .007) 

Valid Real Constants (8 storage locations) : 

21.98753829457168 

1.0000000 

7.9D3 

7.9D03 f (i.e., 7.9 x 10^ = 7900.0) 

7.9D+03 

7.9D+3 

7.9D-03 

7.9D0 

0.0 

7D3 
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(i.e., 7,9 


X 10-3 = .0079) 


(i.e., 7.9 


X 10° = 7.9) 


(i.e., 0.0 


X 10° = 0.0) 


(i.e., 7 X 


103 = 7000) 



Invalid Real Constants: 



3,471.1 

l.E 



7.9D 

1.2E+113 
21.3D90 
23.5E+97 



(missing a decimal point) 

(embedded comma) 

(missing a one- or two-digit integer 

constant following the E. Note that it is not 

interpreted as 1.0 x 10°) 

(missing a one- or two-digit integer 

constant following the D) 

(E is followed by a 3 digit 

integer constant) 

(value exceeds the magnitude permitted; 

that is, 21.3 x 10»o>1663) 

(value exceeds the magnitude permitted; 

that is, 23.5 x 1097>i663) 



COMPLEX CONSTANT 



i Definition 
I- 



Complex Constant - an ordered pair of signed or unsigned real 
constants separated by a comma and enclosed in parentheses. A 
complex constant may assume one of two forms : 

1. From one through seven decimal digits optionally followed by an 
E decimal exponent. In this form, each number in the pair 
occupies four storage locations. 

2. Either one through seven decimal digits followed by a D decimal 
exponent or 8 through 16 decimal digits optionall y followed by a 
D decimal exponent. In this form each number in the pair 
occupies eight storage locations. 

Magnitude: (either form) or 16"^^ through 16 ^^ (i.e., approxi- 
mately ICS) for each real constant in the pair. 



The real constants in a complex constant may be positive, zero, or 
negative (if unsigned, they are assumed to be positive), but they must 
be in the given range. The first real constant in a complex constant 
represents the real part of the complex number; the second represents 
the imaginary part of the complex number. 

Examples : 

Valid Complex Constants: 



(3.2,- 1.86) 

(-5.0E+03,.16E+02) 

(U,0E+03,.16E+02) 

(2.1,0.0) 

(4.70+2,1.9736148) 



(has the value 3.2-1,861) 

(has the value -5000. +16. Oi) 

(has the value 4000. +16. Oi) 

(has the value 2.1+O.Oi) 

(has the value 470. +1. 97361481) 



Where 



i =^^ 
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Invalid Complex Constants : 

(292704,1.6 97) (the real part does not 

contain a decimal point) 

(1. 2E113, 279., 3) (the real part contains 

an invalid decimal exponent) 



LOGICAL CONSTANTS 



r ■ 1 

I Definition | 

|. ^ ^ 

Logical Constan t - There are two logical values : 

.TRUE. 

.FALSE. 

L J 

A logical constant must be preceded and followed by a period. The 
logical constants .TRUE,. and .FALSE. specify that the value of the 
logical variable they replace or the term of the expression they are 
associated with is true or false, respectively. (See the section 
"Logical Expressions.") 



LITERAL CONSTANTS 



r ■ ' 1 

I Definition | 

^ ^ 

I I 

I Literal Constant - a string of alphameric and/or special characters | 
I enclosed in apostrophes. | 

L .^ I J 

The number of characters in the string, including blanks, may not be 
greater than 255. Since apostrophes delimit literal data, a single 
apostrophe within such data is represented by double apostrophes. 

Examples ; 

• DATA* 

•INPUT/OUTPUT AREA NO. 2' 

• X-COORDINATE Y-COORDINATE Z-COORDINATE' 

•3.14» 

'DON'T' 



VARIABLES 

A FORTRAN variable is a symbolic representation of a qucintity that 
may assume different values- The value of a variable may change either 
for different executions of a program or at different stages within the 
program. 

Elements of the Language 13 



For example, in the statement: 

A = 5.0+B 

both A and B are variables. The value of B is determined by some 
previous statement and may change from time to time. The value of A 
varies whenever this computation is performed with a new value for B. 

VARIABLE NAMES 



r •^ 

I Definition | 

^ _^ 

I Variable Name - from 1 through 6 alphameric (i.e., numeric, 0-9, | 
I or alphabetic, A - Z and $) characters, the first of which must be | 
I alphabetic. | 

L .J 

A variable name may not contain special characters (see Appendix A). 
Variable names are symbols used to distinguish one variable from 
another, A name may be used in a source program in one (and only one) 
way (e.g., the name of a variable and that of a subprogram may not be 
identical in the same source program) . 

The use of meaningful variable names can serve as an aid in 
documenting a program. That is, someone other than the programmer may 
look at the program and understand its function. For example, to 
compute the distance a car traveled in a certain amount of time at a 
given rate of speed, the following statement could have been written: 

X = Y * Z 

where * designates multiplication. However, it would be more meaningful 
to someone reading this statement if the programmer had written: 

DIST = RATE + TIME 

Examples : 

Valid Variable Names: 

JOHN 

B292 

VAN 

RATE 

L17N0Y 

SQ704 

Invalid Variable Names : 

B292704 (contains more than six characters) 

UARRAY (first character is not alphabetic) 

SI.X (contains a special character) 



VARIABLE TYPES AND LENGTH SPECIFICATIONS 

The type of a variable corresponds to the type of data the variable 
represents. Thus, an integer variable represents integer data, a real 
variable represents real data, etc. 
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For every type of variable, there is a corresponding standard and 
optional length specification which determines the number of storage 
locations reserved for each variable- The following list shows each 
variable type with its associated standard and optional length: 

Variable Type Standard Optional 

Integer 4 2 

Real U 8 

Complex 8 16 

Logical U 1 

The three ways a programmer may declare the type of a variable are by 
use of the: 

1. Predefined specification contained in the FORTRAN language. 

2. IMPLICIT specification statement. 

3. Explicit speicif ication statements - 

The optional length specification of a variable may be declared only 
by the IMPLICIT or Explicit specification statements. If, in these 
statements, no leingth specification is stated, the standard length is 
assumed (see the section "The Type Statements"). 



TYPE DECLARATION BY THE PREDEFINED SPECIFICATION 

The predefined specification is a convention used to specify vari- 
ables as integer or real, as follows: 

1. If the first character of the variable name is I, J, K, L, M, or N, 
the variable* is integer of standard length. 

2, If the first character of the variable name is any other letter, 
the variable is real of standard length. 

This convention is the traditional FORTRAN method of implicitly 
specifying the type of a variable as being either integer or real. In 
all examples that follow in this publication it is presumed that this 
specification holds, unless otherwise noted. 



TYPE DECLARATION BY THE IMPLICIT SPECIFICATION STATEMENT 

This statement allows a programmer to specify the type of variables 
in much the same way as was specified by the predefined convention. 
That is,, in both, the type is determined by the first character of the 
variable name. However, th$ programmer,, using the IMPLICIT statement, 
has the option of specifying which initial letters designate a particu- 
lar variable type- Further, the IMPLICIT statement is applicable to all 
types of variables — integer, real,, complex,, and logical- 

The IMPLICIT statement overrides the variable type as determined by 
the predefined convention. For example, if the IMPLICIT statement 
specifies that variables beginning with the letters A through M are real 
variables, and variables beginning with the letters N through Y are 
integer variables, then the variable ITEM (which would be treated as an 
integer variable under the predefined convention) is now treated as a 
real variable. Note that variables beginning with the letters Z and $ 
are (by the predefined convention) treated as real variables. The 
IMPLICIT statement is presented in greater detail in the section "Type 
Statements." 
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TYPE DECLARATION BY EXPLICIT SPECIFICATION STATEMENTS 

Explicit specification statements (INTEGER,, REAL,, COMPLEX, and 
LOGICAL) differ from the first two ways of specifying the type of a 
variable, in that an explicit specification statement declares the type 
of a particular variable by its name , rather than as a group of 
variables beginning with a particular character . 

For example, assume: 

1. That an IMPLICIT specification statement overrode the predefined 
convention for variables beginning with the letter I by declaring 
them to be real. 

2. That a subsequent Explicit specification statement declared that 
the variable named ITEM is complex. 

Then, the variable ITEM is complex and all other variables beginning 
with the character I are real. Note that variables beginning with the 
letters J through N are specified as integer by the predefined 
convention. 

The Explicit specification statements are discussed in greater detail 
in the section "Type Statements." 



ARRAYS 

A FORTRAN array is a set of variables identified by a single variable 
name. A particular variable in the array may be referred to by its 
position in the array (e.g., first variable, third variable, seventh 
variable, etc. ) . Consider the array named NEXT, which consists of five 
variables, each currently representing the following values: 

273, 41, 8976,, 59, and 2 

NEXT(l) is the representation of 273 

NEXT (2) is the representation of 41 

NEXT (3) is the representation of 8976 

NEXT (4) is the representation of 59 

NEXT (5) is the representation of 2 

Each variable in this array consists of the name of the array (i.e., 
NEXT) followed by a number enclosed in parentheses, called a subscript. 
The variables which comprise the array are called subscripted variables. 
Therefore, the subscripted variable NEXT(l) has the value 27 3; the 
subscripted variable NEXT (2) HAS THE value 41, etc. 

The subscripted variable NEXT (I) refers to the "Ith" subscripted 
variable in the array, where I is an integer variable that can assume a 
value of 1, 2, 3, 4, or 5. 

To refer to the first element of an array, the array name must be 
subscripted. The array name does not represent the first element. The 
number of subscripts must correspond to the declared dimensionality 
except in the EQUIVALENCE statement- 
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r '■ 1 

I General Form | 

I ^ 

Subscripts - may be one of seven forms: 

V 

c' 

v+c' 

v-c* 

c*v 

c*v+c* 

c*v-c' 

Where: v represents an unsigned, nonsubscripted, integer variable. 

c and c* represent unsigned integer constants. 

L J 

Whatever subscript form is used, its evaluated result must always be 
greater than zero. For example., when reference is made to the 
subscripted variable V(I-2) , the value of I should be greater than 2. 

Examples : 

ARRAY (IHOLD) 
NEXT (19) 
MATRIX (1-5) 
A(5*L) 

W(4*M+3) 

An array may consist of up to seven subscript parameters, separated 
by commas. Thus, the following are valid subscripted variables for 
their corresponding arrays : 

Array Name Sub scripted Variable 

A A(5,, 100, J, K+2) 

TABLE TABLE (1, 1, 1, 1, 1, 1, 1) 

B B(I, J,, K„ L, M, N) 

MATRIX MATRIX(I+2„6*JOB-3,KFRAN) 

Consider the following array named LIST consisting of two subscript 
parameters, the first ranging from 1 through 5, the second from 1 
through 3 : 





Columnl 


(Column2 


Columns 


Rowl 


82 


4 


7 


Row2 


12 


13 


m 


Row3 


91 


1 


31 


Row4 


2U 


16 


10 


Row5 


2 


8 


2 



Suppose it is desired to refer to the number in row 2, column 3; this 
would be: 

LIST (2,3) 

Thus, LIST (2,3) has the value 14 and LIST (U,,l) has the value 24. 

Ordinary mathematical notations might use LIST i,j to represent any 
element, of the array LIST. In FORTRAN, this is written as LIST(I,J), 
where I equals 1,2,3,4, or 5, and J equals 1,2, or 3. 
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As a further example, consider the array named COST,, consisting of 
four subscript parameters. This array might be used to store all the 
premiums for a life insurance applicant, given (1) age, (2) sex,, (3) 
health,, and (4) size of life insurance coverage desired- A code number 
could be developed for each statistic, where lAGE represents age, ISEX 
represents sex,, IHLTH represents health, and ISIZE represents policy 
size desired (see Table 1), 



Table 1. Insurance Premium Codes 



AGE 



SEX 



Age in years 

1-5 
6-10 



96-100 



Code 

IAGE=1 
IAGE=2 



IAGE=20 



HEALTH 



Health 

Poor 
Fair 
Good 
Excellent 



Code 

IHLTH=1 
IHLTH=2 
IHLTH=3 
IHLTH=4 



Sex 

Male 
Fema le 



Code 

ISEX=1 
ISEX=2 



POLICY SIZE 



Dollars 

1,000 

3,000 

5,000 

10,000 

25,000 

50,000 

100,000 



Code 

ISIZE=1 
ISIZE=3 
ISIZE=4 
ISIZE=5 
ISIZE=6 
ISIZE=7 
ISIZE=8 



-H 



Suppose an applicant is 14 years old, male, in good health, and 
desires a policy of $25,000. From Table 1, these statistics can be 
represented by the codes: 



IAGE=3 
ISEX=1 
IHLTH=3 
ISIZE=6 



(11 - 15 years old) 

(male) 

(good health) 

($25,000 policy) 



Thus, COST (3, 1, 3, 6) represents the premium for a policy, given 
the statistics above. Note that "lAGE" can vary from 1 to 20, "ISEX" 
from 1 to 2, "IHLTH" from 1 to 4, and "ISIZE" from 1 to 8. The number 
of subscripted variables in the array COST is the number of combinations 
that can be formed for different ages, sex, health, and policy size 
available - a total of 20x2x4x8 or 1280. Therefore, there may be up to 
1280 different premimums stored in the array named COST. 



DECLARING THE SIZE OF AN ARRAY 



The size of an array is determined by the number of subscript 
parameters of the array and the maximum value of each subscript. This 
information must be given for all arrays before using them in a FORTRAN 
program, so that an appropriate amount of storage may be reserved. 
Declaration of this information is made by a DIMENSION statement,, a 
COMMON statement, or by one of the Explicit specification statements 
(INTEGER, REAL, COMPLEX, and LOGICAL); each is discussed in further 
detail in the section "Specification Statements." 
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ARRANGICMENT OF AlflRAYS IN STORAGE 

Arrays are stored in ascending storage locations, with the value of 
the first of their subscripts increasing most rapidly, and the value of 
the last increasing least rapidly. 

Examples ; 

The array named A, consisting of one subscript parameter which varies 
from 1 to 5,, appears in storage as follows: 

A(l) A(2) A(3) A(4) A(5) 

The array named B,, consisting of two subscript parameters, whose 
first subscript vJiiries over the range from 1 to 5, and second varies 
from 1 to 3, appears in ascending storage locations in the following 
order: 



B(l„l) B(2,l) B(3,l) B(U,1) B(5,l) 



3 



U.B(1,2) B(2,2) B(3»2) B(4/2) B(5,,2)-n 



Ub(1, 



3) B(2,3) B(3,3) B(4,,3) B(5,3) 



Note that B(l,2) and Bd, 3) follow in storage B(5,l) and B(5,, 2), 
respect ively* 

The following list is the order of an array named C, consisting of 
three subscript parameters, whose first subscript varies from 1 to 3, 
second varies from 1 to 2, and third varies from 1 to 3: 

C(l,l,l) C(2,,l,l) C(3,l,,l) C(l,2,l) C(2,2,l) C(3,2,,l)-| 

L^C(1.1,2) C(2,l,2) C(3,l,2) C(l,2,2) C(2,2,2) C(3,2,2)-j 

tc(l,l,,3) C(2,l,3) C(3,l,3) C(l,,2,3) C(2,2,3) C(3,2,3) 

Note that C(l,l,2) and C(l,l,3) follow in storage C(3,,2,,l) and C(3,,2,2),, 
respectively. 



EXPRESS IONS 

Expressions in their simplest form consist of a single constant or 
variable. They may also designate a computation or show a relationship 
between two or more constants and/or variables. Expressions may appear 
in arithmetic and logical assignment statements and in certain control 
statements. 

FORTRAN provides two kinds of expressions: arithmetic and logical. 
The value of an arithmetic expression is always a number whose type is 
integer, real,, or complex. However, the evaluation of a logical 
expression always yields a truth value: .TRUE. or .FALSE.. 



ARITHMETIC EXPRESSIONS 

The simplest arithmetic expression consists of a single constant, 
variable, or subscripted variable of the type integer, real, or complex. 
If the constant, variable, or subscripted variable is of the type 
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integer, the expression is in the integer mode, 
real,, the expression is in the real mode, etc. 



If it is of the type 



Examples ; 

Expression 

3 

I 

3.0 

A 

3,1UD3 



B(2*I) 



(2-0,5.7) 
C 



Type of Quantity 
Integer Constant 
Integer Variable 
Real Constant 
Real Variable 
Real Constant 



Real Variable (Specified 
as such in a Type 
statement) 
Complex Constant 
Complex Variable 
(Specified as such in a 
Type statement) 



Mode of Expression 
Integer 
Integer 
Real 
Real 

Real with eight loca- 
tions of storage re- 
served 

Real with four loca- 
tions of storage re- 
served *♦ 
Complex 
Complex 



In the expression B(2*I), the subscript (2*1), which must always 
represent an integer, does not affect the mode of the expression. That 
is, the mode of the expression is determined solely by the type of 
constant,, variable, or subscripted variable appearing in that expres- 
sion. 

More complicated arithmetic expressions containing two or more 
constants and/or variables may be formed by using arithmetic operators 
that express the computation (s) to be performed. 



Arithmetic Operator s 



The arithmetic operators are as follows: 



Arithmetic Operator 

* 
/ 

+ 



Definition 

Exponentiation 

Multiplication 

Division 

Addition 

Subtraction 



RULES FOR CONSTRUCTING ARITHMETIC EXPRESSIONS ; The following are the 



rules for 
operators : 



constructing arithmetic expressions that contain arithmetic 



1. All desired computations must be specified explicitly. That is, if 
more than one constant, variable, subscripted variable, or subpro- 
gram name (see the section "SUBPROGRAMS") appears in an arithmetic 
expression, they must be separated from one another by an arithme- 
tic operator. For example, the two variables A and B will not be 
multiplied if written: 

AxB or AB or A«B 

If multiplication is desired, then the expression must be written: 

A*B or B*A 

2. No two arithmetic operators may appear in sequence in the same 
expression. For example, the following expressions are invalid: 

A*/B and A*-B 
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However,, in the expression, A*-*B,, if the - is meant to be a minus 
sign rather than the arithmetic operator designating subtraction^ 
then the expression could be written: 

A*(-B) 

In effect^ -B will be evaluated first,, and then A will be 
multiplied with it» (For further uses of parentheses,, see Rule 6.) 

The mode of an arithmetic expression is determined by the type and 
length specification of the variables in the expression. Table 2 
indicates how the mode of variables of different types and lengths 
may be determined using the operators: +„ -, *„ /. 



Table 2, 

,r T 

+ - * / 



Determining the Mode of an Expression Containing 
Variables of Different Types and Lengths 



+- 



INTEGER 
(2) 



INTEGER 
(U) 



REAL 
(4) 



-+- 



REAL 
(8) 



COMPll^X 
(8) 



T 1 

COMPLEX 
(16) 



INTEGER 

(2) 

f + 

INTEGER 
(4) 



Integer 
(2) 

Integer 
(U) 

-I- 



Integer j 
(U) I 

Integer j 
(4) j 
-+- 



Real 

Real 
(4) 



Real 
(8) 

Real 
(8) 



Complex 
(8) 

Complex 
(8) 



Complex 
(16) 
+ ^ 

Complex 
(16) 



REAL 
(4) 



Real 
(4) 



Real 
(4) 



Real 
(4) 



Real 
(8) 



Complex 
(8) 



Complex 
(16) 



REAL 
(8) 



Real 
(8) 



Real 
(8) 



Real 
(8) 



Real 
(8) 



Complex 
(16) 



Complex 
(16) 



COMPLEX 
(8) 



Complex 
(8) 



Complex 
(8) 



Complex 
(8) 



Complex 
(16) 



Complex 
(8) 



Complex 
(16) 



COMPLEX 
(16) 



Complex 
(16) 



Complex i 
(16) I 



Complex 
(16) 



Complex 
(16) 



Complex 
(16) 



Complex 
(16) 



From Table 2 it can be seen that there is a hierarchy of type and 
length specification (see the section "The Type Statements") that 
determines the mode of an expression. For example, complex data 
that has a length specification of 16 when combined with any other 
types of constants and variables results in complex data of length 
16, 

Assume that the type of the following variables has been specified 
as follows: 



Variable Names 
ROOT, E 
A, I, F 
CD 



Type 

Real variable 
Integer variables 
Complex variable 



Length Specificatio n 

4,8 

4,2,2 

16 „ 8 



Then the following examples illustrate how constants and variables 
of differing types may be combined using the arithmetic operators: 
+ * -* /• *: 



Expression 
ROOT* 5 
A+3 

C+2.9D10 
E/F+19 



Mode of Expression 
Real of length 4 
Integer of length 4 
Complex of length 16 
Real of length 8 
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C-18.7E05 



Complex of length 16 
Complex of length 8 



The arithmetic operator denoting exponentiation (i.e.,,**) may only 
be used to combine the types of constants,, variables,, and sub- 
scripted variables shown in Table 3. 



Table 3. Valid Combinations With Respect to the Arithmetic 
Operator ** 

r " ■ 1 

I Base Exponent | 

^ ^ 

(Integer or Real (either length)** Integer or Real (either length)] 

I I 

I Complex (either length) *♦ Integer (either length) | 

L . J 



Assume that the type of the following variables has been specified 
as follows,, and that their length specification is standard: 



Variable Names 
ROOT, E 
A, I, F 
C 



Type 

Real variable 
Integer variables 
Complex variable 



Then the following examples illustrate how constants and variables 
of differing types may be combined using the arithmetic operator,, 
**,. 



Examples : 

Expression 
ROOT** (A+2) 
C**A 
ROOT** I 
I**F 

7.98E21**ROOT 
ROOT** 2, IE 5 
A**E 



Type 

(Real** Integer) 

(Complex** Integer) 

(Real** Integer) 

(Integer** Integer) 

(Real**Real) 

(Real**Real) 

(Integer**Real) 



Result 

(Real) 

(Complex) 

(Real) 

(Integer) 

(Real) 

(Real) 

(Real) 



Order of Computation ; Where parentheses are omitted, or where the 
entire arithmetic expression is enclosed within a single pair of 
parentheses, effectively the order in which the operations are 
performed is as follows: 



Operation 

Evaluation of Functions (see the 

section "Subprograms") 
Exponentiation (**) 

Multiplication and Division (* and /) 
Addition and Subtraction (+ and -) 



Hiera rchy 

1st (highest) 

2nd 
3rd 
Uth 



In addition, if two operators of the same hierarchy (with the 
exception of exponentiation) are used consecutively, the two 
operations are performed from left to right. Thus, the arithmetic 
expression A/B*C is evaluated as if the result of the division of A 
by B was multiplied by C. 
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For example, the expression: 

(A*B/C**I+D) 
is effectivcily evaluated in the order: 

a. C**I Call the result X (exponentiation) 

b. A*B Call the result Y (multiplication) 
c- Y/X Call the result Z (division) 

d, Z+D Final operation (addition) 

For exponentiation the evaluation is from right to left- Thus,, the 
expression: 

is evaluated as follows: 

a. B**C Call the result Z 
b« A**Z Final operation 

Use of P ar entheses ; Parentheses may be used in arithmetic expres- 
sions, as in algebra, to specify the order in which the arithmetic 
operations are to be computed,. Where parentheses are used, the 
expression v«7ithin the parentheses is evaluated before the result is 
used. 

For example,, the expression: 

(B+((A+B)*C)+A**2) 

is effectively evaluated in the order: 

a. (A+B) Call the result X 

b. (X*C) Call the result Y 

c. A** 2 Call the result Z 

d. B+Y+Z Final operations 



LOGICAL EXPRESSIONS 

The simplest form of logical expression consists of a single logical 
constant, logical variable, or logical subscripted variable, the value 
of which is always a truth value (i.e., either ..TRUE. or .FALSE.). 

More complicated logical expressions may be formed by using logical 
and relational operators. These expressions may be in one of the three 
following forms : 

1, Relational operators combined with arithmetic expressions whose 
mode is integer or real. 

2- Logical operators combined with logical constants (.TRUE. and 
.FALSE.), logical variables, or subscripted variables. 

3. Logical operators combined with either or both forms of the logical 
expressions described in items 1 and 2. 

Item 1 is discussed in the following section, "Relational Operators"; 
items 2 and 3 are discussed in the section "Logical Operators." 
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Relational Operators 

The six relational operators, each of which must be preceded and 
followed by a period, are as follows: 

Relational Operator Definition 

.GT. Greater than (» 

.GE. Greater than or equal to (>) 

-LT. Less than (<) 

.LE. Less than or equal to (<) 

,EQ. Equal to (=) 

,NE. Not equal to (*) 

The relational operators express an arithmetic condition which can be 
either true or false. Only arithmetic expressions whose mode is integer 
or real may be combined by relational operators. For example, assume 
that the type of the following variables has been specified as follows : 

Variable Names Type 

ROOT, E Real variables 

A, I, F Integer variables 

L Logical variable 

C Complex variable 

Then the following examples illustrate valid and invalid logical 
expressions using the relational operators. 

Examples; 

Valid Logical Expressions Using Relational Operators: 

(ROOT* A) .GT.E 
A.LT-I 

E**2.7.EQ. (5*R00T+a) 
57.9.LE- (4.7+F) 
.5.GE..9+ROOT 
E.EQ-27.3D+05 

Invalid Logical Expressions Using Relational Operators: 

C.LT.ROOT (Complex quantities may never appear in logical 

expressions) 

C.GE- (2.7,5. 9E3) (Complex quantities may never appear in logical 

expressions) 

L.EQ. (A+F) (Logical quantities may never be joined by 

relational operators) 

E**2,EQ97.1E9 (Missing period immediately after the relational 

operator) 

.GT-9 (Missing arithmetic expression before the rela- 

tional operator) 



Logical Operators 

The three logical operators, each of which must be preceded and 
followed by a period, are as follows. (A and B represent logical 
constants or variables, or expressions containing relational operators.) 
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Logical Operato r 

- NOT.. 

.AND. 
-OR. 



De finition 

.NOT. A - if A is .TRUE., then .NOT. A has the value 
.FALSE.; if A is .FALSE., then .NOT. A has the value 
-TRUE. 

A.AND. B - if A and B are both .TRUE., then A.AND.B 
has the value .TRUE.; if either A or B or both are 
.FALSEi,, then A.AND.B has the value .FALSE, 

A-OR.B - if either A or B or both are .TRUE., then 
A. OR.B has the value .TRUE.; if both A and B are 
.FALSE.,, then A, OR.B has the value .FALfJE. 



Two logical operators may appear in sequence only if the second one 
is the logical operator .NOT.. 

Only those expressions which, when evaluated, have the value .TRUE. 
or .FALSE. may be combined with the logical operators to form logical 
expressions. For example^ assume that the type of the following 
variables has been specified as follows : 



Variable Names Type 

Real variables 
Integer variables 
Logical variables 
Complex variable 

Then the following examples illustrate valid and invalid logical 
expressions using both logical and relational operators. 



ROOT, 


E 


A,, 


Ir 


F 


L, 


w 




C 







Examples ; 

Valid Logical Expressions : 

(ROOT* A, GT- A) , ANDw W 
L, AND,. NOT. (I-GT.F) 

(E+5.9D2,GT.2*E) .OR.L 

. NOT, W, AND, , NOT. L 

L, AND- , NOT- W, OR . I . GT . F 

(A**F,GT, ROOT), AND,, NOT, (I.EQ.E) 

Invalid Logical Expressions: 



A-AND.L 
. OR, W 

NOT. (A-GT.F) 

(C.EQ.D.AND.L 

L.AND. .OR.W 

,AND,L 



(A is not a logical expression) 

(,0R, must be preceded by a logical 

expression) 

(missing period before the logical operator 

. NOT. ) 

(a complex variable may never appear in a 

logical expression) 

(the logical operators .AND. and .OR. must 

always be separated by a logical expression) 

(.AND. must be preceded by a logical 

expression) 



Order of Computation s in Logical Expressions; Where parentheses are 
omitted, or where the entire logical expression is enclosed within a 
single pair of parentheses, the order in which the operations are 
performed is as follows: 
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Operation Hierarchy 

Evaluation of Functions 1st (highest) 

Exponentiation (**) 2nd 

Multiplication and division (* and /) 3rd 

Addition and subtraction (+ and -) 4th 

.LT.,.LE.,-EQ., .NE.,«GT-, .GE. 5th 

-NOT. 6th 

.AND. 7th 

.OR, 8th 

For example,, the expression: 

(A- GT- D**B, AND. . NOT- L. OR, N) 

is effectively evaluated in the following order: 



1. 


D**B 


Call the result 


W 


2- 


A.GT-W 


Call the result 


X 


3. 


.NOT.L 


Call the result 


Y 


4- 


X.AND.Y 


Call the result 


Z 


5, 


Z-OR-N 


Final operation 





(exponentiation) 
(relational operator) 
(highest logical operator) 
(second highest logical operator) 



Use of Parentheses in Logical Expressions: Parentheses may be used in 
logical expressions to specify the order in which the operations are to 
be performed. Where parentheses are used, the expression contained 
within the most deeply nested parentheses (that is,, the innermost pair 
of parentheses) is effectively evaluated first. For example,, the 
logical expression: 

( (I.GT- (B+O) .AND.L) 

is effectively evaluated in the following order: 



1. B+C 

2. I.GT.X 

3. Y, AND.L 



Call the result X 
Call the result Y 
Final operation 



The logical expression to which the logical operator -NOT. applies 
must be enclosed in parentheses if it contains two or more quantities. 
For example, assume that the values of the logical variables A and B are 
.FALSE. and -TRUE-, respectively. Then the following two expressions 
are not equivalent: 

-NOT. (A-OR.B) 
. NOT. A- OR. B 

In the first expression, A.OR.B, is evaluated first. The result is 
.TRUE.; but . NOT, (. TRUE- ) implies .FALSE,. Therefore, the value of the 
first expression is .FALSE.. 

In the second expression, .NOT. A is evaluated first. The result is 
.TRUE.; but .TRUE.,OR-B implies .TRUE., Therefore, the value of the 
second expression is .TRUE,, 
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ARITHMETIC AND LOGICAL ASSIGNMENT STATEMENT 



i General Form j 

^ ^ 

a = b 
Where : 

Notes 



a is any subscripted or nonsubscripted variable. 

b is any arithmetic or logical expression. 

a must be a logical variable if,, and only if^ b is a 
logical expression. 
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The FORTRAN Arithmetic and Logical Assignment statement closely 
resembles a conventional algebraic equation; however, the equal sign of 
the FORTRAN Arithmetic statement specifies replacement rather than 
equivalence. That is,, the expression to the right of the equal sign is 
evaluated, and the resulting value replaces the current value of the 
variable to the left of the equal sign. 

Assume that the type of the following variables has been specified 



as: 

Variable Names 

I, J, W 

A,, B, C, D 

E 

G, H 



Type 

Integer variables 
Real variables 
Complex variable 
Logical variables 



Length Specification 

4,4,2 

4„U,8.8 

8 

U,,4 



Then the following examples illustrate valid arithmetic statements 
using constants, variables,, and subscripted variables of different 
types : 



Statements 
A = B 

W = B 

A = I 

1 = 1 + 1 
E = I++J+D 



A = C*D 



G = -TRUE. 



H = ,NOT,G 



Description 

The value of A is replaced by the current value of B. 

The value of B is truncated to an integer value,, and 
the least significant part replaces the value of W. 

The value of I is converted to a real value, and this 
result replaces the value of A. 

The value of I is replaced by the value of I + 1. 

I is raised to the power J and the result is 
converted to a real value to which the value of D is 
added. This result replaces the real part of the 
complex variable E. The imaginary part of the 
complex variable is set to zero- 

The most significant part of the product of C and D 
replaces the value of A. 

The value of G is replaced by the logical constant 

. TRUE. .. 

If G is .TRUE., the value of H is replaced by the 
logical constant .FALSE.. If G is -FALSE., the value 
of H is replaced by the logical constant .TRUE.. 
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G = 3..GT.I The value of I is converted to a real value; if the 

real constant 3. is greater than this result^, the 
logical constant .TRUE, replaces the value of G. If 
3- is not greater than I, the logical constant 
.FALSE. replaces the value of G. 

E = (1.0,2.0) The value of the complex variable E is replaced by 

the complex constant (1.0,2.0). Note that the state- 
ment E = (A,B) where A and B are real variables is 
invalid. 
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CONTROL STATEMENTS 



Normally, FORTRAN statements are executed sequentially; that is, 
after one statement has been executed, the statement immediately 
following it will be executed. This section discusses the statements 
that may be used to alter and control the normal sequence of execution 
of statements in the program. 



THE GO T O STATEME:NTS 

These statements cause control to be transferred to the statement 
specified by a statement number. There are three GO TO 
statements: Unconditional GO TO, Computed GO TO, and Assigned GO TO- 
Every time the same Unconditional GO TO statement is executed, a 
transfer to the same specified statement is made. However, the Computed 
and Assigned GO TO statements cause control to be transferred to one of 
several statements, depending upon the current value of a particular 
variable. 



Unconditional GO TO Statement 



r • ■ 1 

I General Form | 

^ ^ 

I I 

I GO TO XXX XX I 

I I 

I Where: xxxxx is an executable statement number. | 

L . , J 



This GO TO statement causes control to be transferred to the 
statement specified by the statement number. Every subsequent execution 
of this GO TO statement results in a transfer to that same statement. 

Example ; 

50 GO TO 25 
10 A = B + C 



25 C = E++2 



Explanation: 

In the above example, every time statement numbered 50 is executed, 
control is transferred to the statement numbered 25. 
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Computed GO TO Statement 



i General Form 

^ 

GO TO (Xi, Xa, X3» ...,,Xn)w i 



H 



Where: 



■a.r Xat 



*,,Xn» ciJ^e executable statement numbers. 



i is a nonsubscripted integer variable which is in the 
range: 1 < i < n 

L . _J 



This statement causes control to be transferred to the statement 
numbered x^, Xar X3« . . . , or 3Cn» depending on whether the current value 
of i is 1, 2, 3,-.., or n, respectively. If the value of i is outside 
the allowable range, the next statement is executed. 



Example: 



GO TO (25, 10, 50, 7), ITEM 



50 A = 



25 L = 



B+C 



E**2+A 



C-GT-D.AND.F.LE.G 



10 B = 21.3E02 



Explanation; 



In this example, if the value of the integer variable ITEM is 1, 
statement 2 5 will be executed next. If ITEM is equal to 2, statement 10 
will be executed next, and so on. 
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The ASSIGN and As signed GO TO Statements 



r 1 

I General Form | 

^ „ ^ 

ASSIGN i TO m 



GO TO m, (Xi»X2f X3# • « • f Xn) 

Where: i is an executable statement number. 

Xj., Xa * X3 , . . . .g Xn are executable statement numbers . 

m is a nonsubscripted integer variable of length U to which 
is assigned one of the following statement numbers: 
Xi » X2 f X3 » - « « Xn • 



The Assigned GO TO statement causes control to be transferred to the 
statement numbered Xi^Xa^Xa, . . . ,or Xn» depending on whether the current 
assignment of m is Xj.,X2#X3»- - •#03^ Xn# respectively. For example, in 
the following statement: 

GO TO N, (10, 25, 8) 

If the current assignment of the integer variable N is statement number 
8, the statement numbered 8 is executed next. If the current assignment 
of N is statement number 10, the statement numbered 10 is executed next. 
If N is assigned statement number 25, the statement numbered 25 is 
executed next. 

The current assignment of the integer variable m is determined by the 
last ASSIGN statement executed. Only an ASSIGN statement maiy be used to 
initialize or change the value of the integer variable m. The value of 
the integer m is not the integer statement number; ASSIGN M TO I is not 
the same as I=M. 

Example 1 ; 



ASSIGN 50 TO NUMBER 
10 GO TO NUMBER, (35, 50,, 25,, 12,, 18) 



50 A = B + C 



In the above example,, statement 50 is executed immediately after 
statement 10, 
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Example 2 ; 



ASSIGN 10 TO ITEM 



13 GO TO ITEM, (8, 12, 25, 50, 10) 



8 A = B + C 



10 B = C + D 

ASSIGN 25 TO ITEM 
GO TO 13 



25 C = E**2 



Explanation ; 

In the above example, the first time statement 13 is executed, 
control is transferred to statement 10. On the second execution of 
statement 13, control is transferred to statement 25. 



ADDITIONAL CONTROL STATEMENTS 



Arithmetic IF Statement 



r ~i 

I General Form | 

^ ^ 

IF (a) Xi»2^2»X3 

Where: a is an arithmetic expression which is not complex. 
2£i»2^2»X3 are statement numbers. 

L J 

This statement causes control to be transferred to the statement 
numbered Xi,X2,X3 when the value of the arithmetic expression (a) is 
less than, equal to, or greater than zero, respectively. The first 
executable statement following the arithmetic IF statement should have a 
statement number; otherwise, it can never be referred to or executed. 
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Example : 



IF (A(J,,K)**3-B)10, U, 30 



U D = B + C 



30 C = 



10 E = 



D**2 



(F*B)/D+1 



Explanation ; 

In the above example,, if the value of the expression (A.( J,K) ♦♦3-B) is 

negative, the statement numbered 10 is executed next. If the value of 

the expression is zero, the statement numbered 4 is executed next. If 

the value of the expression is positive, the statement numbered 30 is 
executed next. 



Logica l IF Statemen t 



I General Form | 

^ _ ^ 

IF(a)s 

Where: a is any logical expression. 

s is any statement except a specification stcitement, DO 
statement, or another logical IF statement. 

L . . J 

The logical IF statement is used to evaluate the logical expression 
(a) and to execute or skip statement s, depending on whether the value 
of the expression is .TRUE, or .FALSE., respectively. 

Example 1 : 



5 IF(A.LE-O.O) GO TO 25 

10 C = D + E 

15 IF(A.EQ.B) ANSWER = 2.0+A/C 

20 F = G/H 



25 W = X**Z 
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Explanation : 

In statement 5, if the value of the expression is .TRUE. (i.e., h is 
less than or equal to 0.0), the statement GO TO 25 is executed next, and 
control is passed to the statement numbered 25. If the value of the 
expression is .FALSE. (i.e., A is greater than 0.0), the statement GO 
TO 25 is ignored, and control is passed to the statement numbered 10. 

In statement 15, if the value of the expression is -TRUE. (i.e., A 
is equal to B) , the value of ANSWER is replaced by the value of the 
expression (2,0*A/C), and the statement numbered 20 is executed. If the 
value of the expression is .FALSE. (i.e., A is not equal to B) , the 
value of ANSWER remains unchanged, and the statement numbered 20 is 
executed next. 

Example 2 ; 

Assume that P and Q are logical variables. 



5 IF(P.OR. .NOT.Q)A=B 
10 C = B**2 



Explanation : 

In statement 5, if the value of the expression is .TRUE., the value 
of A is replaced by the value of B and statement 10 is executed next- 
If the value of the expression is .FALSE., the statement A = B is 
skipped and statement 10 is executed. 



DO Statement 



I General Form 
I- 



End of 
Range 



DO 
Variable 



DO 
Where: 



Initial 
Value 



Test 
Value 



Increment 



m 



1113 



X is an executable statement number, that is not defined 
before the DO statement. 

i is a nonsubscripted integer variable. 

5Qa.» IB2(» Esf ^^^ either unsigned integer constants greater 
than zero or unsigned nonsubscripted integer variables whose 
value is greater than zero. The sum m^+ma+l must not exceed 
the size of virtual storage. (ma, is optional; if it is 
omitted, its value is assumed to be 1. In this case, the 
preceding comma must also be omitted.) 



The DO statement is a command to execute repeatedly the statements 
that follow, up to and including the statement numbered x. The first 
time the statements in the range of the DO are executed, i is 
initialized to the value mi; each succeeding time i is increased by the 
value ma. When, at the end of the iteration, i is equal to the highest 



3H 



value that does not exceed ma* control passes to the statement following 
the statement numbered x. Thus, the number of times the statements in 
the range of the DO is executed is given by the expression: 



r 1 

I IL'2 - IDi I 

I I +1 

I ma I 

L J 



where the brackets represent the largest integral value not exceeding 
the value of the expression. If ma is less than m^r the statements in 
the range of the DO are executed once- Upon completion of the DO, the 
DO variable is undefined. 

There are several ways in which looping (repetitively executing the 
same statements) may be accomplished when using the FORTRAN language. 
For example, assume that a manufacturer carries 1000 different machine 
parts in stock,. Periodically, he may find it necessary to compute the 
amount of each different part presently available. This amount may be 
calculated by subtracting the number of each item used, OUT (I), from the 
previous stock on hand, STOCK(I). 

Example; 



5 1=0 

10 1=1+1 

25 STOCK(I)=STOCK(I)- OUT ( I) 

15 IF(I-IOOO) 10,30,30 

30 A=B+C 



Explanation : 

The three statements (5, 10, and 15) required to control the loop 
could be replaced by a single DO statement, as shown in Example 1. 

Example 1 ; 



DO 25 I = 1,1000 
25 STOCK(I) = STOCK(I)-OUT(I) 
3 A=B+C 



Exp l anation : 

In the above example, the DO variable I is set to the initial v alue 

of 1. Before the second execution of statement 25, I is increased by 
the incr e ment 1 and statement 25 is again executed. After 1000 
executions of the DO loop, I equals 1000. Since I is nov7 equal to the 
highest value that does not exceed the test value 10 00, control passes 
out of the DO loop, and statement 30 is executed next. Note that the DO 
variable I is now undefined; its value is not necessarily 1000 or 1001. 
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Example 2 ; 



DO 25 1=1, 10, 2 
15 J=I+K 

25 ARRAY (J) = BRAY (J) 
30 A=B+C 



Explanation ; 



In the above example, statement 25 is the end of the ra nge of the DO 
loop. The DO variable I is set to the initial value of 1. Before the 
second execution of the DO loop, I is increased by the increment 2, and 
statements 15 and 25 are executed a second time. After the fifth 
execution of the DO loop, I equals 9. Since I is now equal to the 
highest value that does not exceed the test value 10, control passes out 
of the DO loop, and statement 30 is executed next. Note that the DO 
variable I is now undefined; its value is not necessarily 9 or 11. 

Programming C onsider ations in Using a DO Loop 

1. The indexing parameters of a DO statment (i, m^ , m^, ma) may not be 
changed by a statement within the range of the DO loop. 



2. 



Range of 
Outer DO 



There may be other DO statements within the range of a DO 
statement. All statements in the range of the inner DO must be in 
the range of the outer DO. A set of DO statements satisfying this 
rule is called a nest of DOS. 

Example 1 ; 

DO 50 I = 1, 4 

A(I) = B(I)**2 

DO 50 J=l, 5 ) 

50 c(j+i) = A(i) ; 

Example 2 ; 

DO 10 INDEX = L, M 

N = INDEX + K 

DO 15 J = 1, 100, 2 
15 TABLE(J) = SUM(J,N)-1 
10 B(N) = A(N) 



Range of 
Inner DO 



Range of 
Outer DO 



Range of 
Inner DO 



3. 



A transfer out of the range of any DO loop is permissible at any 
time. 

If, and only if, a transfer is made from the range of an innermost 
DO loop, transfer back into the range of that innermost DO loop is 
allowed provided none of the indexing parameters Ci»ni3L,m2 /ma) are 
changed outside the range of the DO. A transfer back into the 
range of any other DO in the nest of DOs is not permitted. 
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Example ; 



DO 

DO 







DO 
DO 







■►2 



-:) 



-^ 



5, 



6, 



7, 



Explanation ; 

In the preceding example, the transfers specified by the numbers 
If 2, and 3 are permissible, whereas those specified by 4, 5, and 6 
are not. 

The indexing parameters (i, ,551, ma ,103) may be changed by statements 
outside the range of the DO statement only if no transfer is made 
back into the range of the DO statement using those parameters. 

The last statement in the range of a DO loop (statement x) ^^Y not 
be a GO TO,, Arithmetic IF, PAUSE, STOP, RETURN or another DO 
statement. In addition, the last statement may not be a logical IF 
statement containing any of those statements. 



The use of,, and return from, a subprogram from within any DO 
in a nest of DOs is permitted. 



loop 



CON TINUE Statement 

r ' ' T 

j General Form | 

|. ^ 

I I 

j CONTINUE I 

L J 

CONTINUE is a dummy statement which may be placed anywhere in the 
source program without affecting the sequence of execution. It may be 
used as the last statement in the range of a DO in order to avoid ending 
the DO loop with a GO TO, Arithmetic IF or another DO statement. 

Example 1 ; 



DO 30 I == 1, 20 
7 IF (A(I)-B{I)) 5,30,30 
5 ACI) =A(I) +1.0 

B(I) = B(I) -2.0 

GO TO 7 
30 CONTINUE 
40 C = A(3) + B(7) 
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Explanation ; 

In the preceding example, the CONTINUE statement is used as the last 
statement in the range of the DO to avoid ending the DO loop with the 
statement GO TO 7. 

Example 2 ; 



DO 30 1=1,20 

IF(A(I)-B(I))5,40,40 
5 A(I) = C(I) 

GO TO 3 
HO ACI) = 0.0 
30 CONTINUE 



In Example 2, the CONTINUE statement provides a branch point enabling 
the programmer to bypass the execution of statement 40. 



PAUSE Statement 



r 1 

I General Form | 

^ ^ 

PAUSE 

PAUSE n 

PAUSE ' message ' 

Where: n is an unsigned 1-through 5-digit integer constant. 

messag e is any literal constant. 

L J 

The PAUSE statement causes the program to display 'PAUSE'. If n is 
specified, 'PAUSE n' is displayed; likewise, if ' messag e* is specified, 
'PAUSE message ' is displayed. The program waits until operator inter- 
vention causes it to resume execution, starting with the next statement 
after the PAUSE statement. 



STOP Statement 



r 1 

j General Form j 

j. ^ 

STOP 
STOP n 

Where: n is an unsigned 1-through 5-digit integer constant. 

L J 

This statement terminates the execution of the object program and 
displays n if specified. 



END Statement 



r 1 

I General Form | 

^ ^ 

I I 

I END I 

L : • ^ J 

The END statement is a nonexecutable statement that defines the end 
of a source program or source subprogram for the compiler. Physically, 
it must be the Icist statement of each program or subprogram. 

The END statement must be contained on a single line; however,, 
interspersed blanks between the characters E, N, and D are permitted. 
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INPUT/OUTPUT STATEMENTS 



The input/output statements enable a user to transfer data, belonging 
to a named collection of data,, between I/O devices (such as disk units, 
card readers, and magnetic tape units) and internal storage- The named 
collection of data is called a data set and is not restricted to device 
correspondence. A data set is referred to by an unsigned integer 
constant or integer variable. Formerly, this reference was called a 
symbolic unit number. However, since it more appropriately refers to 
the data rather than any specific device, this number is referred to in 
this publication as the data set reference number. 

For the FORTRAN user, a data set is considered to be a continuous 
string of data which may be subdivided into FORTRAN records. This 
subdivision of data sets into FORTRAN records is stated by the use of 
one or more of the following: 

1. A FORMAT statement referred to by an I/O statement 

2. An I/O list appearing in an I/O statement 

3. A NAMELIST name appearing in an I/O statement 

In addition to subdividing data sets into records, a FORMAT statement 
may be used to declare the form in which the data is to be transmitted. 

There are five I/O statements: READ, WRITE, END FILE, REWIND, and 
BACKSPACE. The READ and WRITE statements are used to transfer data into 
or from internal storage. The END FILE statement defines the end of a 
data set; the REWIND and BACKSPACE statements control the positioning of 
data sets. 

Even though the I/O statements described below are device indepen- 
dent, in that a given I/O statement may be applicable to a data set on 
any number of devices or device types , it is often meaningful to 
consider the original source, or ultimate destination of the data being 
transferred- Thus, for the sake of demonstration, subsequent examples 
will be in terms of card input and print- line output. 



BASIC IN PUT/OUTPUT STATEMENTS 

The basic input/output statements are READ and WRITE- The statements 
FORMAT and NAMELIST, though not I/O statements, may be used in 
conjunction with certain forms of READ and WRITE statements. All four 
statements are presented in greater detail in the following sections. 
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READ STATEMENT 



J. . ^ 

I General Form | 

^ ^ _ ^ 

READ (a, b„ END=c„ ERR=d) list 

Where: a is an unsigned integer constant or an integer variable of 
length 4 that represents a data set reference number. 

b is either the statement number or array name of the FORMAT 
statement describing the data being read, or a NAMELIST 
name. 

c is the statement number to which transfer is made upon 
encountering the end of the data set. 

d is the statement number to which transfer is made upon 
encountering an error condition in data transfer., 

list is a series of variable or array names, separated by 
commas, which may be indexed and incremented. They specify 
the number of items to be read and the locations in storage 
into which the data is placed. 
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The READ Statement may take many different forms. For example, the 
parameters END=c and ERR=d are optional and, therefore, may or may not 
appear in a READ statement. Furthermore, either the lis t or the 
parameter b may be omitted. 

When one or more of the parameters END=c or ERR=d are used after the 
a and b portion of a READ statement, they may appear in any order within 
the parentheses. For example, the following are valid READ statements: 

READ(5,50,ERR=10)A,B,C 
READ(5, 25, END=15) D,E,F,F,H 
READ (N, 30, ERR=100, END=8) X,,Y,Z 

If a transfer is made to a statement specified by the END parameter, 
no indication is given the program as to the number of items in the list 
(if any) read before encountering the end of the data set. If an END 
parameter is not specified in a READ statement, the end of the data set 
terminates execution of the object program. 

If a transfer is made to a statement specified by the ERR parameter, 
no data is read into the list items associated with the record in error. 
No indication is given the program as to which input record or records 
are in error; only that one or more data items read into the list may be 
in error,. If an ERR parameter is not specified in a READ statement, an 
error terminates execution of the object program. 

The three basic forms of the READ statement are: 

READ (a,x) 
READ ( a, b) list 
READ (a) list 

The parameters END=c and ERR=d may be used in the combination 
described above in each of these three forms. 
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The Form READ (a,, x) 

This form is used to read data from the data set associated with a 
into the locations in storage specified by the NAMELIST name x. The 
NAMELIST name x is a single variable name that refers to a specific list 
of variables or array names into which the data is placed. A specific 
list of variable or array names receives a NAMELIST name by use of a 
NAMELIST statement. The programmer need only use the NAMELIST name in 
the READ (a,x) statement to reference that list thereafter in the 
program. 

The format and rules for constructing and using the NAMELIST 
statement are described in the following text. 

J. . . ^ 

I General Form | 

^ ^ 

NAMELIST/x/a,br .. .•c/y/d,e, rf/z/2,h, .. . , i 

Where: x#y# and z, . . . are NAMELIST names. 

a,b,c,d, ... are variable or array names. 

L ; . J 

The following rules apply to defining and using a NAMELIST name: 

1. A NAMELIST name consists of from 1 through 6 alphameric characters, 
the first of which is alphabetic. 

2. A NAMELIST name is enclosed in slashes. The list of variable or 
array names belonging to a NAMELIST name ends with a new NAMELIST 
name enclosed in slashes or with the end of the NAMELIST statement. 

3. A variable name or an array name may belong to one or more NAMELIST 
names . 

U. A NAMELIST name may be defined only once by its appearance in a 
NAMELIST statement and must be so defined before its use. After it 
is defined in the NAMELIST statement, the NAMELIST name may appear 
only in input or output statements thereafter in the program. 

5. A NAMELIST statement may appear anywhere in a FORTRAN program prior 
to its use in a READ/WRITE statement. 

6. Variable or array names appearing anywhere in a NAMELIST statement 
or NAMELIST name may not appear in a FUNCTION, SUBROUTINE, or ENTRY 
statement. 

Example : 

Assume that A, I, and L are array names. 



NAMELIST /NAM1/A,B,I,J,L/NAM2/A,C,J,K 



READ (5,NAM1) 
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Explanation ; 

The above READ statement causes the record that contains the input 
data for the variables and arrays that belong to the N/^ELIST name 
referenced, NAMl, to be read from the data set associated v^ith the data 
set reference number 5. 

When a READ statement references a NAMELIST name, input data in the 
form described in the following text is read from the designated input 
data s et . 

Input Data 

The first character in the record is always ignored. The second 
character of the first record of a group of data records to be read must 
be g, immediately followed by the NAMELIST name. This name is followed 
by any combination of data items 1 and 2 below, separated by commas. (A 
comma after the last item is optional.) 

The form the data items may take is : 

1. Var iable nam e = constant 

The variable name may be a subscripted variable name or a single 
variable name. Subscripts must be integer constants. 

2- Array name = set of constants (separated by commas) 

The set of constants may be in the form "k* constant," where k is 
an unsigned integer called the repeat constant. It represents the 
number of successive elements in the array to be initialized by the 
specified constant. The number of constants must be equal to the 
number of elements in the array. 

Constants used in the data items may be integer, real, literal, 
complex, or logical data. If the constants are logical data, they may 
be in the form T or .TRUE,. and F or .FALSE.. 

Any selected set of variable or array names belonging to the NAMELIST 
name appearing on the first record may be used as specified by items 1 
and 2 in the preceding text. Names that are made equivalent to these 
names may not be used unless they also belong to the NAMELIST name. 

The end of a group of data is signaled by the character string &END 
with no embedded blanks and all appearing in the same record. 

Blanks must not be embedded in a constant or repeat constant, but may 
be used freely elsewhere in a data record. The last item on each record 
that contains data items must be a constant followed by a comma. (The 
comma is optional on the record that precedes the record containing 
SEND.) 

Example ; 

Assume that L is an array consisting of one subscript parameter 
ranging from 1 to 10. 
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Column 2 

t 
First Data Card: 6NAM1 
Second Data Card: I (2, 3) =5, J=4, 
Third Data Card: A(3)=4.0, L=2,3,8*4, 



Last Data Card: SEND 

Explanation : 

If this data is input to be used with the NAMELIST and READ 
statements previously illustrated, the following actions take place. 
The first data card is read and examined to verify that its name (and 
the data items that follow) is consistent with the NAMELIST name in the 
READ statement. If that NAMELIST name is not found it reads to the next 
namelist group. When the second data card is read, the integer 
constants 5 and 4 are placed in 1(2,3) and J, respectively. When the 
third data card is read, the real constant U.O is placed in A(3) . Also, 
since L is an array not followed by a subscript, the entire array is 
filled with the succeeding constants. Therefore, the integer constants 
2 and 3 is placed in L(l) and L(2), respectively, and the integer 
constant 4 is placed in L(3) , L( 4) , . . . , L(10) . 



The Form READ (a, b) List 

This form is used to read data from the data set associated with a 
into the locations in storage specified by the variable names in the 
list. The list , used in conjunction with the specified FORMAT statement 
b (see the section "FORMAT statement"), determines the number of items 
(data) to be read, the locations, and the form the data will take in 
storage. 

Example 1 : 

Assume that the variables A, B, and C have been declared as integer 
variables. 



75 FORMAT (GlO, G8, G9) 



READ (J, 75) A, B, C 



Explanation : 

The above READ statement causes input data from the data set 
associated with data set reference number J to be read into the 
locations A, B, and C according to the FORMAT statement referenced 
(statement 75). That is, the first 10 positions of the record are read 
into storage location A; the next 8 positions are read into storage 
location B; and the next 9 positions are read into storage location C. 

The list can be omitted from the READ (a,b) list statement. In this 
case, a record is skipped or data is read from the data set associated 
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with a into the locations in storage occupied by the FORMAT statement 
numbered b. 

Example 2 : 



98 FORMAT ("HEADING') 



READ (5,98) 



Explanation ; 

The above statements would cause the characters H,, E,, A, D„ I, N, and 
G in storage to be replaced by the next 7 characters in the data set 
associated with data set reference number 5, 

Example 3 ; 



98 FORMAT (GIO- • HEADING* ) 



READ (5,98) 



Explanation ; 

The above statements would cause the next record in the data set 
associated with data set reference number 5 to be skipped. No data is 
transferred into internal storage because there is no list item that 
corresponds with format code GIO. 

The Form READ (a) Lis t 

The form READ (a) lis t of the READ statement causes binary data 
(internal form) to be read from the data set associated with a into the 
locations of storage specified by the variable names in the list. Since 
the input data is always in internal form, a FORMAT statement is not 
required. This statement is used to retrieve the data v/ritten by a 
WRITE (a) list statement- 

Example; 

READ (5) A, B,, C 

Explanation; 

This statement causes the binary data from the data set associated 
with data set reference number 5 to be read into the storage locations 
specified by the variable names A„ B,, and C. 

The l ist may be omitted from the READ (a) list statement. In this 
case, a record is skipped. 
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Example ; 

READ (5) 

Explanat-ion ; 

The above statements would cause the next record in the data set 
associated with data set reference number 5 to be skipped. No data is 
transferred into internal storage. 



Indexing I/O Lists 

Variables within an I/O list may be indexed and incremented in the 
same manner as those within a DO statement. These variables and their 
indexes must be included in parentheses. For example,, suppose it is 
desired to read data into the first five positions of the array A. This 
may be accomplished by using an indexed list as follows; 

15 FORMAT (G10.3) 



READ (2,15) (A(I), 1=1,5) 
This is equivalent to: 
15 FORMAT (G10,3) 



DO 12 I = 1,5 
12 READ (2,15) A(I) 

As with DO statements, a third indexing parameter may be used to 
specify the amount by which the index is to be incremented at each 
iteration. Thus, 

READ (2,15) (A(I), 1=1,10,2) 

causes transmission of values for A(l), A(3), A(5), A(7),, and A(9). 

Furthermore, this notation may be nested. For example, the state- 
ment: 

READ (2,15) ( (C(I,J),,D(I,J),J=1,3),,I=1„4) 

would transmit data in the following order: 

C(l.l), D(l,l),, C(l,2). D(l,2), C(l,3), D(l.,3) 

C(2,l), D(2^1),, C(2,2), D(2,,2), C(2,3), D(2,3) 

C(3,l),, D(3,l), C(3,,2), D(3,2), CC3,,3), D(3,3) 

C(4,l), D(U,1), C(U.2). D(4,2), C(4,,3), D(U,3) 

Since J is the innermost index, it varies more rapidly than I. 

As another example, consider the following: 

READ (2,25) I, (C ( J) , J=l,, I) 

The variable I is read first and its value then serves as an index to 
specify the number of data items to be read into the array C. 
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If it is desired to read data into an entire arrciy, it is not 
necessEiry to index that array in the I/O list. For example, assume that 
the array A consists of one subscript parameter varying in the range of 
1 to 10. Then the following READ statement referring to FORMAT 
statement numbered 5 : 

READ (2,5) A 

would cause data to be read into A(l)# A(2),...,A(10). 

The indexing of I/O lists applies to WRITE lists, as well as READ 
lists. 



Reading Format Statements 

FORTRAN provides the facility for variable FORMAT statements by 
allowing a FORMAT statement to be read into an array in storage and 
using the data in the array as the FORMAT specifications for subsequent 
I/O statements. 

For example, the following statements result in A„ B,, and the array C 
being read, converted, and stored according to the FORMAT specifications 
read into the array FMT at object time: 

DIMENSION FMT (18) 
1 FORMAT (18AU) 
READ (5,1) FMT 
READ (5, FMT) A,B,, (C ( I) , 1=1, 5) 

1. The name of the variable FORMAT specification must appear in a 
DIMENSION statement, even if the array size is only 1. 

2. The form of the format codes read into the FMT array at object time 
must take the same form as a source program FORMAT statement,, 
except that the word FORMAT is omitted (see the section "The FORMAT 
Statement") . 



WRITE STATEMENT 



r ■ 1 

I General Form | 

I ^ 

WRITE (a, b) list 

Where: a is an unsigned integer constant or an integer variable of 
length U that represents a data set reference number. 

b is either the statement number or array name of the FORMAT 
statement describing the data being written, or a NAMELIST 
name. 

list is a series of variable or array names, separated by 

commas, which may be indexed and incremented. They specify 

the number of items to be written and the locations in 
storage from which the data is taken. 
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The WRITE statement may take many different forms. For example, the 
list or the parameter b may be omitted. 

The three basic forms of the WRITE statement are: 

WRITE (a, x) 
WRITE (a,, b) list 
WRITE (a) list 



The Form WRITE (a, x ) 

This form is used to write data from the locations in storage 
specified by the NAMELIST name x into the data set associated with a 
(see the section "The Form READ(a,x)"). 

Example ; 

WRITE(6,NAM1) 

Explanation ; 

This statement causes all variable and array names (as well as their 
values) that belong to the NAMELIST name, NAMl, to be written on the 
data set associated with data set reference number 6. 

When a WRITE statement references a NAMELIST name: 

Ik, All variables and arrays and their values belonging to the NAMELIST 
name will be written out, each according to its type. The complete 
array is written out by columns. 

2, The output data will be written such that; 

a. The fields for the data will be large enough to contain all the 
significant digits. 

b. The output can be read by an input statement referencing the 
NAMELIST name. 

Example : 

Assume that A is a 3 by 3 array. 



NAMELIST/NAMl/A, B, I , D 
WRITE (8,NAM1) 



Assuming that the output is punched on cards, the format would be: 

Column 2 

t 
First Output Card; SNAMl 

Second Output Card: A=3,4, 4.5, 6.2, 25.1, 
Third Output Card: 9.0, -15.2,-7.6, 0.576Ebl2, 

Fourth Output Card: 2, 717, B=3. 14,1=10,0=0. 378E-15, 
Fifth Output Card: SEND 
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The Form WRITE (ci,b) List 

This form is used to write data in the data set associated with a 
from the locations in storage specified by the variable neimes in the 
list . The list, used in conjunction with the specified FORMAT statement 
b, determines the number of items (data) to be written^ the locations, 
and the form the data will take in the data set. 

Example? 1 ; 

In the following example, assume that the variables A, B, and C have 
been declared as integer variables. 

75 FORMAT (C510, G8„ G9) 



WRITE (J, 75) A, B, C 

Explanation ; 

The above WRITE statement causes output data to be written in the 
data set associated with the data set reference number J, from the 
locations A, B,, C, according to the FORMAT statement referred to 
(statement 75). That is,, the 10 rightmost digits in A are written in 
the data set associated with the data set reference number J; the next 8 
positions in the data set will contain the 8 rightmost digits in B; and 
the next 9 positions in the data set will contain the 9 rightmost digits 
in C. 

The list may be omitted from the WRITE (a,fb) list statememt. In this 

case, a blank record is inserted, or data is written in the data set 

associated with a from the locations in storage occupied by the FORMAT 
statement b- 

Example 2 : 

98 FORMAT C ' HEADING*) 



WRITE (5,98) 

The above statements would cause a blank and the characters H, E, A, 
D, I, N, and G in storage to be written in the data set associated with 
data set reference number 5. 

Example 3 ; 

98 FORMAT (GlO, 'HEADING') 



WRITE (5,98) 

Explanation ; 

The above statements would cause a blank record to be placed in the 
data set associated with data set reference number 5 . No data is 
transferred into the data set. 
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The Form WRITE (a) List 



The WRITE (a) list form of the WRITE statement causes binary data 
(internal form) from the locations of storage specified by the variable 
names in the list to be written in the data set associated with a. 
Since the output data is always in internal form, a FORMAT statement is 
not required. The READ (a) list statement is used to retrieve the data 



written by a WRITE (a) list statement. 



Example: 



WRITE (5)A, B, C 



Explanation; 

The statement causes the binary data from the locations specified by 
the variable names A, B, and C to be written in the data set associated 
with data set reference number 5. 



FORMAT STATEMENT 



I General Form 
I- 



xxxxx FORMAT ( c^. , Ca, . . . , Cn/Ci ' » Ca ' » . . . # Cn '/...) 

Where: xxxxx is a statement number (1 through 5 digits ) - 

Ci,C2#--.fCn and Ci* , Ca* r • • • r Cn' are format codes which may 
be delimited by one of the separators: comma, slash, or 
parenthesis. These codes specify the length, decimal point 
(if any), and position of the data in the data set. 



/ may be used to separate FORTRAN records. 



The FORMAT statement is used in conjunction with the READ and WRITE 
statements in order to specify the desired form of the data to be 
transmitted. The form of the data is varied by the use of different 
format codes. The twelve format codes are: G, T, X, P, literal. A, I, 
F, E, D, H, and L. Any number used in a FORMAT statement except a 
statement number or a literal must be less than 256. 

USE OF THE FORMAT STATEMENT: This section contains general information 
on the FORMAT statement. The points discussed below are illustrated by 
the examples that follow. 

1, FORMAT statements are nonexecutable and may be placed anywhere in 
the source program, 

2, A FORMAT Statement may be used to define a FORTRAN record, as 
follows: 

a. If no slashes or additional parentheses appear within a FORMAT 
statement, a FORTRAN record is defined by the beginning of the 
FORMAT statement (left parenthesis) to the end of the FORMAT 
statement (right parenthesis). Thus, a new record is read when 
the format control is initiated (left parenthesis); a new 
record is written when the format control is terminated (right 
parenthesis) , 
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Example !: 



L corresponds to 1 

FORTRAN record 

b. If slashes appear within a FORMAT statement, FORTRAN records 
are defined by the beginning of the FORMAT statement to the 
first slash in the FORMAT statement, from one slash to the next 
succeeding slash, or from the last slash to the end of the 
FORMAT statement. Thus, a new record is read whem the format 
control is initiated,, and thereafter a record is read upon 
encountering a slash; a new record is written upon encountering 
a slash or when format control is terminated. 



Example ; 

XX XX X FORMAT ( / / ) 

< > < > < > 

I I I 



L each corresponds to 

1 FORTRAN record 

c. If more than one level of parentheses appears within a FORMAT 
statement, a record is defined by the beginning of the FORMAT 
statement to the end of the FORMAT statement; thereafter, from 
the first- level left parenthesis from the right of the FORMAT 
statement to the end of the FORMAT statement. 

Example 1 ; 

1 2 21 
xxxxx FORMAT ( ( ( )) ) 



•each corresponds to 
1 FORTRAN record 



Example 2 ; 



Oil 110 
xxxxx FORMAT ( ( ) ( ) ) 



> 



< > 

I 



-each corresponds to 
1 FORTRAN record 
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When defining a FORTRAN record by a FORMAT statement, it is 
important to consider the original source (input) or ultimate 
destination (output) of the record. For example,, if a FORTRAN 
record is to be punched for output, the record should not be 
greater than 8 characters. For input, the FORMAT statement should 
not define a FORTRAN record longer than the record referred to in 
the data set, 

3, Blank output records may be introduced or input records may be 
skipped by using consecutive slashes (/) in a FORMAT statement. If 
there are n consecutive slashes at the beginning or end of a FORMAT 
statement, n input records are skipped or n blank records are 
inserted between output records, respectively. If n consecutive 
slashes appear anywhere else in a FORMAT statement, the number of 
records skipped or blank records inserted is n-1. 

4, Successive items in an I/O list are transmitted according to 
successive format codes in the FORMAT statement, until all items in 
the list are transmitted. If there are more items in the list than 
there are codes in the FORMAT statement, control transfers to the 
preceding left parenthesis of the FORMAT statement and the same 
format codes are used again with the next record. If there are 
fewer items in the list, the remaining format codes are not used. 

5, A format code may be repeated as many times as desired by preceding 
the format code with an unsigned integer constant. 

6, A limited parenthetical expression is permitted to enable repeti- 
tion of data fields according to certain format codes within a 
longer FORMAT statement. Two levels of parentheses,, in addition to 
the parentheses required by the FORMAT statement,, are permitted. 
The second level of parentheses facilitates the transmission of 
complex quantities. • 

7, When transferring data on input or output, the type of format code 
used, type of data, and type of variables in the I/O list should 
correspond. 

8- In the following examples, the output is shown as a printed line. 
A carriage control character 'x' (see "Carriage Control") is 
specified in the FORMAT statement but does not appear in the first 
print position of the print line. This carriage control character 
appears as the first character of the output record on any I/O 
medium except the printed line. 

G Format Code 



r 

I General Form 

j. 

aGw.s 

Where: a is optional and is an unsigned integer constant used to 
denote the number of times the same format code is repeti- 
tively referenced. 

w is an unsigned integer constant less than or equal to 255 
specifying the total field length. 

s is an unsigned integer constant specifying the number of 
significant digits. 
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The G format code is a generalized code, in that it may be used to 
determine the desired form of data, whether it be integer, real, 
complex, or logical. 

The .s portion may be omitted when transferring integer or logical 
data. If present, it is ignored. When real data is transferred, the w 
portion of the G format code includes four positions for a decimal 
exponent field. 

If the real data, say n, is in the range 0.1<n<10**s, where s is the 
s portion of the format code Gw.s^, this exponent field is blank. 
Otherwise, the r€;al data is transferred with an E or D decimal exponent, 
depending on the length specification (either 4 or 8 storage locations, 
respectively) of the real data. 

For simplification, the following examples deal with the printed 
line; however, the concepts developed apply to all input/output media. 

Example 1 ; 

Assume that the variables A, B, C, and D are of type real, v^hose values 
are 292,7041, 82.43441, 136,7632, .8081945, respectively, 

1 FORMAT C •X*,G12,4,G12,5,G12,4,G12.7) 

2 FORMAT C'xV,Gl3,4,Gl3.5,Gl3,4) 

3 FORMAT (•x',Gl3,4) 



WRITE (5, b) A, B, C, D 



Explanation : 

a. If b had been specified as 1, the printed output would be: (b 
represents a blank) 

Print Position 1 Print Position 4 8 

t t 

bbb292,7bbbbbb82.434bbbbbbbl36.7bbbb.8081945bbbb 

b. If b had been specified as 2, the printed output would be: 

Print Position 1 Print Position 39 

t t 

bbbb292,7bbbbbbb8 2.434bbbbbbbbl36.7bbbb Line 1 

bbb0,8081bbbb Line 2 

From the example, it can be seen that by increasing the field 
width reserved (w) , blanks are inserted. 

c. If b had been specified as 3, the printed output would be: 

Print Position 1 

t 

bbbb292,7bbbb Line 1 

bbbb82.43bbbb Line 2 

bbbbl36,7bbbb Line 3 

bbbO.SOBlbbbb Line 4 
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From the example, it can be seen that the same format code was 
used for each variable in the list . Each repetition of the 
same format code caused a new line to be printed. 

Example 2 : 

Assume that the variables I, J, K, and L are of type integer, whose 
values are 292, 443U28, 4908081, and 40018, respectively. 

1 FORMAT ( 'x',GlO,2G7,G5) 

2 FORMAT CxV.Ge) 

3 FORMAT (•x»,4G10) 



WRITE (5, b) I, J, K, L 

Explanation : 

a. If b had been specified as 1, the printed output would be: 

Print Position 1 Print Position 29 

t t 

bbbbbbb292b443U28490808140018 Line 1 

The same results would be achieved, if FORMAT statement 1 had 
been written: 

FORMAT (•x',GlO, G7, G7, G5) 

Note that the .s portion of the G format may be omitted when 
transmitting integer data. 

b. If b had been specified as 2, the printed output would be: 

Print Position 1 

f 

bbb292 Line 1 

U43U28 Line 2 

908081 Line 3 

bU0018 Line 4 

Note that the second format code G6 is an incorrect specifi- 
cation for the third variable K, i.e., 4908081. Thus, the 
leftmost digit is lost. In general, when the width specifi- 
cation w is insufficient, the leftmost characters are not 
printed. 

c. If b had been specified as 3, the printed output would be: 

Print Position 1 Print Position 40 

t t 

bbbbbbb292bbbb443428bbb4908081bbbbb40018 Line 1 

From the above example, it can be seen that increasing the 
field width w improves readability. 
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Example 3 ; 

Assume that the variable I is integer (length 2),, A and B are real 
(length 4), D is real (length 8), C is complex (length 8), and L is 
logical (length 1) whose values are 292, 471.93, 81.91, 6.9310072, 
(2.1,3.7), and ,TRUE.» respectively. 

1 FORMAT ("x\G3,2G9.2,Gl3,7,2G8,2,G3) 

2 FORMAT ("x\G3/*x',2Gl0.2/'x',G9.1/'x',2G8.2,G3) 

3 FORMAT (//•x*,G3,2G9.2//'x' ,G13.7,2G8.2,G3///) 



WRITE (5,b) I,A,B,D,C,L 

Explanation ; 

a. If b has been specified as 1, the printed output vrauld be: 

Print Position 1 Print Position 53 

t t 

292b0.47Eb03bb81ibbbbb6.931007bbbbb2.1bbbbb3.7bbbbbbT 

When complex data is being transmitted, two format codes are 

required. The real and imaginary parts are each treated as 

separate real numbers, and the parentheses and comma are not 
printed as part of the output. 

b. If b has been specified as 2, the printed output v/ould be: 

Print Position 1 

292 Line 1 

bb0.47Eb03bbb81.bbbb Line 2 

bbb6bbbb Line 3 

b2.1bbbbb3.7bbbbbbT Line 4 

From the example, it can be seen that the use of the slash (/) 
to separate two format codes causes the data, not yet printed, 
to be printed on a new line. If the output data is to be 
punched on cards, the slash specifies that the following data 
will be punched on another card. 

c. If b has been specified as 3, the printed output would be: 



Print Position 1 

t 

(blank line) 

(blank line) 

292b0 . 47Eb03bb81. bbbb 

(blank line) 

b6,931007bbbbb2.1bbbbb3.7bbbbbbT 

(blank line) 

(blank line) 

(blank line) 



Line 1 
Line 2 
Line 3 
Line 4 
Line 5 
Line 6 
Line 7 
Line 8 



In the example, note that 2 consecutive slashes appearing at 
the beginning and 3 at the end of the series of format codes 
causes blank lines to be inserted as shown^ However, the two 
consecutive slashes appearing elsewhere in the FORMAT state- 
ment cause the insertion of a blank line, as shown in line 4. 
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The principles illustrated in the previous output examples also apply 
when using the READ statement on input. In addition, there are further 
considerations when using the FORMAT statement on input or output. 

1. When reading real input data with a G format code, a decimal point 
must be included. 

2. The use of additional parentheses (up to two levels) within a 
FORMAT statement is permitted to enable the user to repeat the same 
format code when transmitting data. For example, the statement: 

10 FORMAT (2(G10,6,G7.1) ,G4) 

is equivalent to: 

10 FORMAT (GIO.^, G7 . 1 , G10.6, G7.1, G4) 

3. If the data exists with a D decimal exponent, it is transferred 
with this D decimal exponent. 

4. If a multiline listing is desired such that the first two lines are 
to be printed according to a special format and all remaining lines 
according to another format, the last format code in the statement 
should be enclosed in a second pair of parentheses. For example, 
in the statement: 

FORMAT (•x*,G2,2G3.1/'x* ,Gl0.8/("x' ,3G5.1) ) 

If more data items are to be transmitted after the format codes 
have been completely used, the format repeats from the last left 
parenthesis. Thus, the printed output would take the form: 

G2,G3.1,G3.1 
G10.8 

G5.1,G5.1,G5.1 
G5.1,G5.1,G5.1 



As another example, consider the statement: 

FORMAT ( •x*,G2/2('x' , G3 ,G6 . 1) ,G9 . 7) 

If there are thirteen data items to be transmitted, the printed 
output on a WRITE statement would take the form: 

G2 

G3,G6.1, •xV,G3,G6.1,G9.7 

G3,G6.1, 'x' ,G3,G6.1,G9.7 

G3,G6.1 



Numeric Format Codes (I,F ,E,D) 

Four types of format codes are available for the transfer of numeric 
data. These are specified in the following form: 
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r ' 1 

I General Form | 

|. _ ^ 

aFw,d 
aEw . d 
aDw.d 

Where: a is optional and is an unsigned integer constant used to 
denote the number of times the same format code is repeti- 
tively referenced. 

I, F, E, and D are format codes- 

w is an unsigned integer constant less than or equal to 255 
specifying the total field length of the data. 

d is an unsigned integer constant specifying the number of 
decimal places to the right of the decimal point, i.e., the 
fractional portion. 



For purposes of simplification, the following discussion of format 
codes deals with the printed line. The concepts developed apply to all 
input/output media. 



I Format Code 

The I format code is used in conjunction with the trcinsferral of 
integer data. The code 110 may be used to print integer delta; 10 print 
positions are reserved for the number. It is printed in this 
10-position field right- justified (that is, the units position is at the 
extreme right). 

If the number to be transmitted is greater than 10 positions, the 
excess leftmost digits are lost. If the number has less than 10 digits, 
the leftmost print positions are filled with blanks. If the quantity is 
negative, the position preceding the leftmost digit contains a minus 
sign. In this case,, an additional position should be specified in the 
total field length for the minus sign. On input, if the field length 
specification w is greater than the number of digits being jread into a 
field, the integer data is right- justified and high-order zeros are 
inserted. 

The following examples show how each of the quantities on the left is 
printed according to the format code 13 (b represents a blank) : 

Internal Value Printed Value 

721 721 

-721 721 (incorrect because of insufficient 

specification) 
-12 -12 

568im 114 (incorrect because of insufficient 

specification) 
bbO 

-5 b-5 

9 bb9 
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F Format Code 

The F format code is used in conjunction with the transferral of real 
data that does not contain a decimal exponent. For F format codes,, w is 
the total field length reserved, and d is the number of places to the 
right of the decimal point (the fractional portion) . This differs from 
the G format code» where the number of significant digits is specified. 
The total field length reserved must include sufficient positions for a 
sign (if any) and a decimal point. The sign,, if negative,, is printed. 

If insufficient positions are reserved by d, the fractional portion 
is truncated from the right. If excessive positions are reserved by d, 
zeros are filled in on the right. The integer portion of the number is 
handled in the same fashion as numbers transmitted by the I format code. 

The following examples show how each of the quantities on the left is 
printed according to the format code F5.2: 

Internal Value Printed Value 

12.17 12,17 

-41-16 41,16 (incorrect because of insufficient 

specification) 
-,2 -0,20 

7.3542 b7.35 (last two digits of accuracy lost 

because of insufficient specification) 
-1. -1.00 

9,03 b9.03 

187,64 87,64 (incorrect because of insufficient 

specification) 



D and E Format Code s 

The D and E format codes are used in conjunction with the transferral 
of real data that contains a D or E decimal exponent,, respectively. A D 
format code indicates a field length of 8; an E format code indicates a 
field length of 4. For D and E format codes, the fractional portion is 
again indicated by d. The w includes field d, spaces for a sign, the 
decimal point, and four spaces for the exponent. (For output,, space 
should be reserved for at least one digit preceding the decimal point.) 

The exponent is the power of 10 by which the number must be 

multiplied to obtain its true value. The exponent is written with a D 

or an E, followed by a space for the sign and two spaces for the 
exponent (maximum is 75) , 

The following examples show how each of the quantities on the left is 
printed according to the format codes (DIO, 3/ElO, 3) : 

Internal Value Printed Value 

238. b0.238Db03 

-.002 -0.200E-02 

.00000000004 bO,400D-10 

-21.0057 -0.210Eb02 (Last three digits of accuracy 

lost because of insufficient 

specification) 

When reading input data, the start of the exponent field must be 
marked by an E, or, if that is omitted, by a + or - sign (not a blank) . 
Thus, E2„ E+2, +2, +02, E02, and E+02 all have the same effect and are 
permissible decimal exponents for input. 
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Numbers for E, D,, and F format codes need not have their decimal 
point punched. If it is not present, the decimal point is supplied by 
the d portion of the format code. If it is present in the card,, its 
position overrides the position indicated by the d portion of the format 
code. 



L Format Code 



I General Form 
|. __. 



aLw 



Where; 



a is optional and is an unsigned integer constant used to 
denote the number of times the same format code is 
repetitively referenced- 

w is an unsigned integer constant less than or equal to 255, 
specifying the number of characters of data. 



Lw. 



Logical variables may be read or written by means of the format code 



On input, the first T or F encountered in the next w characters of 
the input record causes a value of .TRUE. or .FALSE., respesctively, to 
be assigned to the corresponding logical variable. If the field w 
consists entirely of blanks,, a value of .FALSE. is assumed. 

On output, a T or an F is inserted in the output record corresponding 
to the value of the logical variable in the I/O list. The single 
character is preceded by w - 1 blanks. 



A Format Code 



i General Form 
k 



aAw 



Where: 



a is optional and is an unsigned integer constant used to 
denote the number of times the same format code is repeti- 
tively referenced. 

w is an unsigned integer constant less than or equal to 255,, 
specifying the number of characters of data. 



The format code Aw is used to read or write data. If w is equal to 
the number of characters corresponding to the length specification of 
each item in the I/O list, w characters are read or written., 

On input, if w is less than the length specification of each item in 
the I/O list, w characters are read and the remaining right-most 
characters in the item are replaced with blanks. If w is greater than 
the length specification, the number of characters equal to the 
difference between w and the length specification are skipped and the 
remaining characters are read,. 
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On output, if w is less than the length specification of the item in 
the I/O list, the printed line will consist of the leftmost w characters 
of the item. If w is greater than the length specification, the printed 
line will consist of the characters right- justified in the field and 
will be preceded by blanks. Therefore it is important to always 
allocate enough area in storage to handle the characters being written 
(see the section "The Type Statements"). 

Exa m ple 1 : 

Assume that the array ALPHA consists of one subscript parameter ranging 
from 1 through 20, The following statements could be written to "copy" 
a record from one data set to another whose ultimate destination is a 
card punch. 



10 FORMAT (20A4) 



READ (5,10) (ALPHA( I) ,1=1,20) 



WRITE (6,10) (ALPHA (I), ,1=1, 20) 



Explanation ; 

In example, the READ statement would cause 20 groups of characters to 
be read from the data set associated with data set reference number 5. 
Each group of four characters would be placed into the 20 positions in 
storage starting with ALPHA (1) and ending with ALPHA (20). The WRITE 
statement would cause the 20 groups of four characters to be written on 
the data set associated with data set reference number 6. 

Example 2 ; 

As another example, consider all the variable names in the list of 
the following READ statement to have been explicitly specified as REAL, 
and the array CONST to have been specified as having one subscript 
parameter ranging from 1 through 10. Then assuming the following input 
data is associated with data set reference number 5: 

ABCDE XYZ$1234567890b 

where . . . represents the alphabetic characters F through W and b means 
a blank, the following statements could be written: 



10 FORMAT (27A1,10A1,A1) 
20 FORMAT ( 'x ' , 6 (7A1, 5X) ) 



READ (5,10)A,B,C,D,E,F,G,H,I, 

1 J,K,L,M,N,0,P,Q,R,. 

2 S,T,U,V,W,X,Y,Z,$. 

3 (CONST (IND),IND=1, 10), BLANK 
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DO 50 INDEX =1,5 



WRITE (6,20) G, R, O, U, P, BLANK, CONST ( INDEX) , 

1 B, L, O, C, K, BLANK , CONST ( INDEX) ,, 

2 F, I, E,L,,D, BLANK, CONST (INDEX),, 

3 G, R,, O, U , P , BLANK, CONST ( INDEX+5 ) ,, 
U B, L, O, C , K, BLANK, CONST ( INDEX+5 ) , 
5 F, I, E, L, D, BLANK, CONST ( INDEX+5 ) 



50 CONTINUE 



Explanation ; 

The READ statement would cause the 37 alphameric characters and the 
blank in the data set associated with data set reference number 5 to be 
placed into the storage locations specified by the variable names in the 
READ list. Thus, the variables A through Z receive the values A through 
Z, respectively; the variable $ receives the value $; the numbers 1 
through 9, and 0, are placed in the 10 fields in storage starting with 
CONST(l) and ending with CONST(IO); and the variable BLANK receives a 
blank. The WRITE statement within the DO loop would cause the following 
heading to be printed. A subsequent WRITE statement within the DO loop 
could then be written to print the corresponding output data. 



Print Position 1 
t 



Print Position 67 
t 



GROUP 1 BLOCK 1 FIELD 1 GROUP 6 BLOCK 6 FIELD 6 

(output data) 
GROUP 2 BLOCK 2 FIELD 2 GROUP 7 BLOCK 7 FIELD 7 

(output data) 



GROUP 5 BLOCK 5 FIELD 5 GROUP 

(output data) 



BLOCK 



FIELD 



Literal Data in a Format Statement 



Literal data consists of a string of alphameric and special charac- 
ters written within the FORMAT statement and enclosed in apostrophes. 
The string of characters must be less than or equal to 255. For 
example: 

25 FORMAT (' 1964 INVENTORY REPORT') 

An apostrophe character within the string is represented by two 
successive apostrophes. For example, the characters DON'T are rep- 
resented as: 
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DON* '• T 

The effect of the literal format code depends on whether it is used 
with an input or output statement. 

INPUT 

A number of characters, equal to the number of characters between the 
apostrophes, are read from the designated data set. These characters 
replace, in storage., the characters within the apostrophes. 

For example, the statements: 



FORMAT ( • HEADINGS ' ) 



READ (3,5) 



would cause the next nine characters to be read from the data set 
associated with data set reference number 3; these characters would 
replace the blank and the eight characters H,E,A,Dw I,N,G, and S in 
storage. 

OUTPUT 

All characters (including blanks) within the apostrophes are written 
as part of the output data. Thus, the statements: 



FORMAT ( • THIS IS ALPHAMERIC DATA* ) 



WRITE (2,5) 



would cause the following record to be written on the data set 
associated with the data set reference number 2: 

THIS IS ALPHAMERIC DATA 



H Format Code 



r 1 

I General Form | 

1^ ^ 

wH 

Where: w is an unsigned integer constant less than or equal to 255, 
specifying the number of characters following H. 

L J 
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The H format code is used in conjunction with the transferral pf 
literal data. 

The format code wH is followed in the FORMAT statement by w (w<25$) 
characters,. For example, 

5 FORMAT (31H THIS IS ALPHAMERIC INFORMATION) 

Blanks are significant and must be included as part of the count w- 
The effect of wH depends on whether it is used with input or output. 

1, On input, w characters are extracted from the input record aiid 
replace the w characters of the literal data in the FORMAT 
statement. 

2, On output, the w characters following the format code are writtfen 
as part of the output record. 



X Format Code 



r ' ^n 

I General Form ; | 

^ . ^ ^^ 

wX 

Where: w is an unsigned integer constant less than or equal to 255, 
specifying the number of blanks to be inserted on output or 
the number of characters to be skipped on input. 

L , , ^J 

When the wX (w<255) format code is used with a READ statement (i.eL, 
on input), w characters are skipped before the data is read in. FJDr 
example,, if a card has six 10-column fields of integer quantities, ahd 
it is not desired to read the second quantity, then the statement: 

5 FORMAT (I10„10X, 4110) 

may be used, along with the appropriate READ statement. 

When the wX format code is used with a WRITE statement (i.e., pn 
output) , w characters are left blank. Thus, the facility for spacing 
within a printed line is available. For example, the statement: 

10 FORMAT ( •xr,3(F6,2,5X)) 

may be used with an appropriate WRITE statement to print a line as 
follows : 

123. U5bbbbb817, 32bbbbb52a, 67bbbbb 
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T Format Code 



r 

I General Form 

F 



Tw 

Where: w is an unsigned integer constant less than or equal to 255, 
specifying the position in a FORTRAN record where the 
transfer of data is to begin . 

Input and output may begin at any position by using the format code 
Tw (w<255) . Only when the output is printed does the correspondence 
between w and the actual print position differ. In this case,, because 
of the. carriage control character, the print position corresponds to 
w-1, as in the following example: 

5 FORMAT (TUO, '1964 INVENTORY REPORT' T80,, 'DECEMBER' Tl, ' PART 
NO, 10095') 

The above FORMAT statement would result in a printed line as follows: 

Print Print Print 

Position 1 Position 39 Position 79 

t t ♦ 

PART NO. 10095 1964 INVENTORY REPORT DECEMBER 

The statements : 

5 FORMAT (T40, ' HEADINGS') 



READ (3,5) 

would cause the first 39 characters of the input data to be skipped, and 
the next 9 characters would then replace the blank and the characters 
H,E|,A,D,I,N,G and S in storage. 

The T format code may be used in a FORMAT statement with any type of 
format code. For example, the following statement is valid: 

5 FORMAT (TlOO, FlO.3, T50, E9.3, Tl, ' ANSWER IS') 



Scale Factor - P 

The representation of the data, internally or externally, may be 
modified by the use of a scale factor followed by the letter P preceding 
a format code. 

The scale factor is defined for input and output as: 

external quantity = internal quantity x 10 scale factor 

For input, when scale factors are used in a FORMAT statement, they 
have effect only on real data which does not contain an E or D decimal 
exponent. For example, if input data is in the form xx.xxxx and, it is 
desired to use it internally in the form -xxxxxx, the format code used 
to effect this change is 2PF7.4. 
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INPUT 

As another example, consider the input data: 

27bbb-93. 2094bb-175« 8041bbbb55« 3647 : 

where b represents a blank. 
The statements: 

5 FORMAT (I2,3F11.4) 

READ (6,5) K,,A,,B,C 

would cause the variables in the list to assume the values : 

K : 27 B : -175.80U1 

A : -93.2094 C : 55.3647 

The statements : 

5 FORMAT (I2,1P3F11.4) 

READ (6,5) K,A,B,C 

would cause the variables in the list to assume the values: 

K : 27 B : -17.5804 

A : -9.3209 C : 5.5364 

The statements: 

5 FORMAT (I2,-1P3F11.4) 

READ (6,5) K„A,B,C 

would cause the variable in the list to assume the values : 

K : 27 B : -1758. 041x 

A : -932.094X C : 553.647x 

where the x represents an extraneous digit, 

OUTPUT 

Assume that the variables K,A,B, and C have the values: 

K : 27 B : -175.8041 

A : -93.2094 C : 55.3647 

the statements: 

5 FORMAT (I2,1P3F11.4) 

WRITE (4,5) K,A,B,C 
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would cause the variables in the list to output the values: 

K : 27 B : -1758. 041x 

A : -932.094X C : 553.647x 

where the x represents an extraneous digit. 

The statements: 

5 FORMAT (I2,-1P3F11.4) 



WRITE (4,5) K,A,B,C 

would cause the variables in the list to output the values: 

K : 27 B : -17.5804 

A : -9.3209 C : 5.5364 

For output, when scale factors are used, they have effect only on real 
data. However, this real data may contain an E or D decimal exponent. 
A positive scale factor used with real data which contains an E or D 
decimal exponent increases the number and decreases the exponent. Thus, 
if the real data was in a form using an E decimal exponent, and the 
statement FORMAT (1X,I2, 3E13. 3) used with an appropriate WRITE statement 
resulted in the printed line: 

27bbb-0.932Eb02bbb-0.175Eb03bbbb0.553Eb02 

Then the statement FORMAT (1X,I2,1P3E13. 3) used with the same WRITE 
statement would result in the printed output: 

27bbb-9.320Eb01bbb-1.75 8Eb02bbbb5.536Eb01 

The scale factor is assumed to be zero if no other value has been 
given. However, once a value has been given, it will hold for all 
format codes (i.e., those that correspond to real data) following the 
scale factor within the same FORMAT statement. This also applies to 
format codes enclosed within an additional pair of parentheses,. Once 
the scale factor has been given, a subsequent scale factor of zero in 
the same FORMAT statement must be specified by OP. 



Carriage Control 

When records written under format control are prepared for printing,, 
the following convention for carriage control applies: 

First Character Carriage Advance Before Printing 

Blank One Line 

Two lines 

1 To first line of the next page 
+ No advance 

The first character of the output record may be used for carriage 
control and is not printed. It appears in all other media as data. 
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ADDITIONAL INPU T/OUTPOT STATEMENTS 

The Statements END FILE, REWIND, and BACKSPACE are used to control 
the data sets, as described in the following text. 



END FILE Statement 



r " • 1 

I General Form | 

|. ^ 

END FILE a 

Where: a is an unsigned integer constant or integer variable of 

length 4 that represents a data set reference numbejr. 

L . . J 

The END FILE statement defines the end of the data set associated 
with a. A subsequent WRITE statement defines the beginning of a new 
data set. 



REWIND Statement 



r 1 

I General Form j 

|. _ _ ^ 

REWIND a 

Where: a is an unsigned integer constant or integer Vciriable of 
length 4 that represents a data set reference numbesr. 

L , J 

The REWIND statement causes a subsequent READ or WRITE statement 
referring to a to read data from or write data into the first data set 
associated with a.. 



BAq:KSPACE Statemen t 



r — " T 

I General Form | 

^_^ ^ ^ 

Backspace a 

Where: a is an unsigned integer constant or integer variable of 
length U that represents a data set reference number, 

L J 

The BACKSPACE statement causes the data set associated with a to 
backspace one record. If the data set associated with a is already at 
its beginning, execution of this statement has no effect. 
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SPECIFICATION STATEMENTS 



The specification statements provide the compiler with information 
about the nature of the data used in the source program. In addition, 
they supply the information required to allocate locations in storage 
for this data. Specification statements describing data may appear 
anywhere in the source program, but must precede any statements which 
refer to that data. 



THE TYPE STATEMENTS 

There are two kinds of type statements: the IMPLICIT specification 
statement and the Explicit specification statements (INTEGER, REAL, 
COMPLEX, and LOGICAL). 

The IMPLICIT specification statement enables the user to: 

1. Specify the type of a group of variables or arrays according to the 
initial character of their names. 

2. Specify the amount of storage to be allocated for each variable 
according to the associated type. 

The Explicit specification statements enable the user to: 

1- Specify the type of a variable or array according to their 
particular name. 

2- Specify the amount of storage to be allocated for each variable 
according to the associated type. 

3. Specify the dimensions of an array. 

4- Assign initial data values for variables and arrays- 



IMPLICIT Statement 



r —1 

I General Form | 

^ _ __^ 

IMPLICIT type*s (a±,a2, ...),... , type *s ia±,a2, . . . ) 

Where: type represents one of the following: INTEGER, REAL, COM- 
PLEX, or LOGICAL. 

*s is optional and represents one of the permissible length 
specifications for its associated type. 

§.11 3.2,..- represent single alphabetic characters each 
separated by commas, or a range of characters (in alphabetic 
sequence) denoted by the first and last characters of the 
range separated by a minus sign (e.g., (A-D) ) . 

L . J 
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The IMPLICIT statement', if specified, should be the first statement 
in a main program, and the second statement in a FUNCTION,, SUBROUTINE,, 
or BLOCK DATA subprogram. 

The IMPLICIT type statement enables the user to declare the type of 
the variables appeiaring in his program (i.e., integer, real, complex, or 
logical) by specifying that variables beginning with certain designated 
letters are of a certain type. Furthermore, the IMPLICIT statement 
allows the programmer to declare the number of locations to be allocated 
for each in the group of specified variables. The type a variable can 
assume, along with the permissible length specifications are as follows: 

Type Length Specification 

INTEGER 2 or U (standard length is 4) 

REAL 4 or « (standard length is 4) 

COMPLEX 8 or 16 (standard length is 8) 

LOGICAL' 1 or 4 (standard length is 4) 

For each type there is a corresponding standard length speicif ication. 
If this standard length specification (for its associated type) is 
desired, the *s may be omitted in the IMPLICIT statement. That is, the 
variables will assume the standard length specification. For each type 
there is also a corresponding optional length specification., If this 
optional length specification is desired,, the *s^ must be included within 
the IMPLICIT statement. 

Example 1 ; 

IMPLICIT REAL (A-H,, 0-Z,$)„ INTEGER (I-N) 

Explanation ; 

All variables beginning with the characters I through N are declared 
as INTEGER. Since no length specification was explicitly given (i.e., 
the *s^ was omitted),, four storage locations (the standard length for 
INTEGER) are allocated for each variable. 

All other variables (those beginning with the characters A through H, 
O through Z, and $) are declared as REAL with four storage locations 
allocated for each. 

Note that the statement in Example 1 performs the exact same function 
of typing variables as the predefined convention (see "Type Declaration 
by the Predefined Specification"). 

Example 2 : 

IMPLICIT INTEGER*2(A-H) , REAL*8(I-K), LOGICAL (L, M, N) 

Ex planation ; 

All variables beginning with the characters A through H are declared 
as integer, with two storage locations allocated for each. All 
variables beginning with the characters I through K are declared as 
real, with eight storage locations allocated for each. All variables 
beginning with the characters L, M, and N are declared as logical, with 
four locations allocated for each. 

Since the remaining letters of the alphabet (O through Z and $) were 
left undefined by the IMPLICIT statement, the predefined convention will 
take effect. Thus, all variables beginning with the characters O 
through Z and $ are declared as real, each with a standard length of 
four locations. 
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Example 3 : 

IMPLICIT C0MPLEX*16(C-F) 

Explanation ; 

All variables beginning with the characters C through F are declared 
as complex^ each with eight storage locations reserved for the real part 
of the complex data and eight storage locations reserved for the 
imaginary part. The types of the variables beginning with the charac- 
ters A, B,, G through Z, and $ are determined by the predefined 
convention. 



Explicit Specification Statements 



r 

I General Form 

I- 



Type *s a*Si(ki)/Xa./,b*S2(k2)/X2/# • • • »z*Sn(kn)/Xn/ 
Where: Type is INTEGER, REAL, LOGICAL, or COMPLEX. 

*s» *Si, *S2, ,, . , ♦Sn are optional. Each s represents one of 
the permissible length specifications for its associated 
type . 

a,b»,..^z represent variable, array, or function names (see 
the section, "SUBPROGRAMS") 

(kj.) » ^ISa^ » • • • » ^iSn) ai^e optional. Each k is composed of 1 
through 7 unsigned integer constants, separated by commas, 
representing the maximum value of each subscript in the 
array. Each k may be an unsigned integer variable only when 
it appears in a Type statement in a subprogram. 

/Xi./»/X2/» . . . #/Xn/ are optional and represent initial data 
values. 

The Explicit specification statements declare the typ e (INTEGER, 
REAL, COMPLEX, or LOGICAL) of a particular variable or array by its 
name, rather than by its initial character . This differs from the other 
ways of specifying the type of a variable or array (i.e., the predefined 
convention and the IMPLICIT statement) . In addition, the information 
necessary to allocate storage for arrays (dimension information) may be 
included within the statement. However, if this information does not 
appear in an Explicit specification statement, it must appear in a 
DIMENSION or COMMON statement (see "DIMENSION Statement" or "COMMON 
Statement") . 

Initial data values may be assigned to variables or arrays by use of 
/Xn/» where Xn is a constant or list of constants separated by commas. 
This set of constants may be in the form "r* constant", where r is an 
unsigned integer, called the repeat constant. No element may have more 
than one initial value given in the same program. A function name may 
not have an initial value assigned to it. An initially defined variable 
or a variable of an array may not be in blank common. In a labeled 
common block, they may be initially defined only in a BLOCK DATA 
subprogram. 
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In the same manner in which the IMPLICIT statement overrides the 
predefined convention, the Explicit specification statements override 
the IMPLICIT and predefined convention. If the length specification is 
omitted (i.e.,*s), the standard length per type is assumed. 

Example 1 ; 

INTEGER+2 ITEM/76/,, VALUE 

Explanation ; 

This statement declares that the variables ITEM and VALUE are of type 
integer, each with two storage locations reserved. In addition, the 
variable ITEM is initialized to the value 76. 

Example 2 ; 

COMPLEX C,D/(2-l,U.7)/,E*16 

Explanation ; 

This statement declares that the variables C, D, and E are of type 
complex. Since no length specification was explicitly given, the 
standard length is assumed. Thus, C and D each have eight storage 
locations reserved (four for the real part, four for the imaginary part) 
and D is initialized to the value (2-1,4.7). In addition, 16 storage 
locations are reserved for the variable E. Thus, if a length specifi- 
cation is explicitly written, it overrides the assumed standard length. 

Example 3 ; 

REAL*8 ARRAY, HOLD, VALUE+4, ITEM(5,5) 

Explanation ; 

This statement declares that the variables ARRAY, HOLD, VALUE, and 
the array named ITEM are of type real. In addition, it declares the 
size of the array ITEM. The variables ARRAY and HOLD have eight storage 
locations reserved for each; the variable VALUE has four storage 
locations reserved; and the array named ITEM has 200 storage locations 
reserved (eight for each variable in the array). Note that when the 
length is associated with the type (e.g., REAL*8), the length applies to 
each variable in the statement unless explicitly overridden (as in the 
case of VALUE* 4),. 

Example U ; 

REAL A(5,5)/20*6.9E2,5*1.0/, B (100)/100*0. 0/,TOAD*8 (5) /5*0. 0/ 

Explanation ; 

This statement declares the size of each array, A and B, and their 
type (real). The array A has 100 storage locations reserved (four for 
each variable in the array) and the array B has UOO storage locations 
reserved (four for each variable). In addition, the first 20 variables 
in the array A are initialized to the value 6,9E2 and the last five 
variables are initialized to the value 1.0. All 100 variables in the 
array B are initialized to the value 0.0, The array TOAD has 40 storage 
locations reserved (eight for each variable). In addition, each 
variable is initialized to the value 0.0. 
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Adjustab l e Dimensions 

As shown in the previous examples, the maximum value of each 
subscript in an array was specified by a numeric value- These numeric 
values (maximum value of each subscript) are known as the absolute 
dimensions of an array and may never be changed. However, if an array 
is used in a subprogram (see the section "Subprograms") and is not in 
Common, the size of this array does not have to be explicitly declared 
in the subprogram by a numeric value. That is, the Explicit specifi- 
cation statement, appearing in a subprogram, may contain integer 
variables that specify the size of the array. When the subprogram is 
called, these integer variables then receive their values from the 
calling program. Thus, the dimensions (size) of a dummy array appearing 
in a subprogram are adjustable and may change each time the subprogram 
is called. 

The absolute dimensions of an array must be declared in a calling 
program. The adjustable dimensions of an array, appearing in a 
subprogram, should be less than or equal to the absolute dimensions of 
that array, as declared in the calling program. 

The following example illustrates the use of adjustable dimensions: 

Example ; 

Calling P rogram Subprogram 



REAL* 8 A(5,5) 



SUBROUTINE MAPMY ( . . . , R, L, M, . . . ) 



CALL MAPMY(.. . ,A,2,3,. . .) 



REAL* 8. ..,R(L,M) ,. 



DO 100 1=1, L 



Explanation; 

The statement REAL*8 A(5,5) appearing in the calling program declares 
the absolute dimensions of the array A. When the subroutine MAPMY is 
called, the dummy argument R assumes the array name A, and the dummy 
arguments L and M assume the values 2 and 3, respectively. The 
subscripted variables of the array A appearing in the calling program 
occupy unique storage locations in the following order: 

A(l,l) A(2,l) A(3,l) A(U,1) A(5,l) 

A(l,2) A(2,2) A(3,2) A(4,2) A(5,2) 

A(l,3) A(2,3) A(3.3) A(U,3) A(5,3) 

A(l,4) A(2,4) A(3,4) A(4,U) A(5,U) 

A(l,5) A(2,5) A(3,5) A(4,5) A(5,5) 

Thus, in the calling program the subscripted variable A (1,2) refers to 
the sixth subscripted variable in the array A. However, in the 
subprogram MAPMY the subscripted variable Ad, 2) refers to the t hir d 
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subscripted varicible in the array A, namely, A(3,l). This is so because 
the dimensions of the array A as declared in the subprogram are not the 
same as those in the calling program. 

If the absolute dimensions in the calling program were the same as 
the adjusted dimensions in the subprogram, the subscripted variables 
A(l„l) through At 5, 5) in the subprogram would always refer to the same 
storage locations as specified by the subscripted variables A(l,l) 
through A(5,5) in the calling program, respectively. 

The numbers 2 and 3, which became the adjusted dimension of the dummy 
array R, could also have been variables in the argument list of the 
calling program. For example, assume that the following statement 
appeared in the calling program: 

CALL MAPMY ( , A, I,, J, . . . ) 

Then as long as the values of I and J were previously determined,, the 
arguments may be variables. In addition, the variable dimension size 
may be passed through more than one level of subprograms. For example, 
within the subprogram MAPMY could have been a call statement to another 
subprogram in which dimension information about A could have been 
passed. 

If any dimension of an array is variable, that dimension and the 
array name must be dummy variables (i.e., they must appear in a 
FUNCTION, SUBROUTINE, or ENTRY statement) - 



ADDITIONAL SPECIFICATION STATEMENTS 



DIMENSION Statement 



r ' • 1 

I General Form | 

■ ^ ^ 

DIMENSION ai(ki) ,a2(k2) , £3 (kg), »an(kn) 

Where: a^, a^.* a3»---# ^n are array names. 

]Si» ]S2» iS3#*--*»lSn are each composed of 1 through 7 unsigned 
intege^r constants, separated by commas, representing the 
maximum value of each subscript in the array. k^ through kn 
may he integer variables of length 4 only when they appear 
in a DIMENSION statement within a subprogram. 

L , . J 

The information necessary to allocate storage for arrays used in the 
source program may be provided by the DIMENSION statement. The 
following examples illustrate how this information may be declared. 



DIMENSION A (10) , ARRAY (5,5,5,5,5), LIST(10,100) 
DIIVIENSION BC25, 50) , TABLE ( 25 , 25, 25) 
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COMMON Statement 



j General Form 
I- 



COMMON /r/a (ka.),b(k2), /r/cCkg) ,d(k ), 

Where: a,b, . . . ,c,d. . , are variable or array names, 

iSi'iSat • • 'iiafk ••• 3^6 optional and are each composed of one 
through seven unsigned integer constants, separated by 
commas, representing the maximum value of each subscript in 
the array, 

/r/... represent optional common block names consisting of 
one through six alphameric characters, the first of which is 
alphabetic. These names must always be embedded in slashes. 



Although the COMMON statement may be used to provide dimension 
information, adjustable dimensions may never be used. 

Variables or arrays that appear in a calling program or a subprogram 
may be made to share the same storage locations with variables or arrays 
in other subprograms by use of the COMMON statement. For example, if 
one program contains the statement: 

COMMON TABLE 
and a second program contains the statement: 

COMMON LIST 

the variable names TABLE and LIST refer to the same storage locations. 

If the main program contains the statements: 

REAL A, B,C 
COMMON A,B,C 

and a subprogram contains the statements : 

REAL X,,Y*Z 
COMMON X,Y, Z 

A shares the same storage location as X; B shares the same storage 
location as Y; and C shares the same storage location as Z. 

Consider the following examples: 

Example 1 ; 



Calling Program 



Subprogram 

SUBROUTINE MAPMY (...) 



COMMON A, B, C, R(IOO) 
REAL A, B,C 
INTEGER R 



COMMON X,, Y, Z, S(IOO) 
REAL X,Y,Z 
INTEGER S 



CALL MAPMY ( ) 
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Ex planation t 



In the calling program, the statement COMMON A,B,C,R(100) would cause 
412 storage locations (four locations per variable) to be reserved in 
the following order: 



Beginning j 
of common | 
area | 

1" 



4 locations 



B 

U locations 



4 locations 



Layout of 
storage 



R(l) 

U locations 



R(IOO) 

U locations 



The statement COMMON X, Y^ Z,, S(IOO) would then cause the variables 
X, Y, Z, and S(1).-,S(100) to share the same storage space as A, B, C, 
and R(1).,.R(100), respectively. 

From the above example, it can be seen that COMMON statements may be 
used to serve an important function: namely, as a medium to implicitly 
transmit data from the calling program to the subprogram. That is, 
values for X„ Y» Z, and S(l) . . .S (100) , because they occupy the same 
storage locations as A, B,, C, and R(l) « , .R(IOO) , do not have to be 
transmitted in the argument list of a CALL statement. 

Example 2 : 

Assume COMMON is defined in a main program and three subprograms as 
follows : 



Main program 
Subprogram 1 
Subprogram 2 
Subprogram 3 



COMMON 
COMMON 
COMMON 
COMMON 



A,B,C 
D,E,F 
Q^R,S,T,U 
V^W.X,Y.Z 



Further,, assume the length specifications of the above variables are so 
defined that the common area is shared as follows: 



~ 


— — . 




T 


~ 


" 


~ 


^^ 


T~ 




— __ 


A 








B 








1 


c 




8 


locations 




^ _4. 


U 


locations 






1 
J. 


2 


locations 








— ^ — ^ 
















D 








E 








1 


F 




8 


locations 

^. 


..«--. 


+- 


H 


locations 

T- 





^ 


1 


2 


locations 


Q 


1 


R 




S 


1 


T 




1 


U 




4 


locations) 
i- 


H 


locations | 
— 4 


2 


locations ] 
_ _ 1- 


2 


locati 


ons| 


2 


locations 




T 




T 




T 






T 






V 


1 


W 




X 


1 


Y 




1 


Z 




U 


locations | 


U 


locations | 


2 


locations | 


2 


locations ] 


2 


locations 





i. 





_ JL_ 




^ X- 








x_ 









In this case, the variables A,B,C and D,E,F may be 
in their respective programs, as may Q, R,S„T 
addition, all programs may validly refer to C,F,U, 
possible to cross-reference D in Subprogram 1 and 
2« Such correspondences are highly data dependent 
may be useful. For instance, if D is a complex 
are real variables, Q and R correspond to the real 
of D, respectively. However, each such cross ref 
mer must be considered on its own merits. 



validly referred to 
,U and V,W„X,Y,Z. In 
and Z- It is also 
Q and R in Subprogram 
and in certain cases 
variable, and Q and R 
and imaginary parts 
erence by the program- 
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Blank and Labeled Common 

In the preceding two examples,, the common storage area (common block) 
established is called a blank common area. That is, no particular name 
was given to that area of storage. The variables that appeared in the 
COMMON statements were assigned locations relative to the beginning of 
this blank common area. However, variables and arrays may be placed in 
separate common areas. Each of these separate areas (or blocks) is 
given a name consisting of one through six alphameric characters (the 
first of which is alphabetic) ; those blocks which have the same name 
occupy the same storage space. 

Those variables that are to be placed in labeled (or named) common 
are preceded by a common block name enclosed in slashes. For example, 
the variables A, B, and C will be placed in the labeled common area HOLD 
by the following statement: 

COMMON/ HOLD/A, B, C 

In a COMMON statement, blank common may be distinguished from labeled 
common by preceding the variables in blank common by two consecutive 
slashes or, if the variables appear at the beginning of the common 
statement, by omitting any block name. For example, in the following 
statement : 

COMMON A, B, C /ITEMS/ X, Y, Z / / D, E, F 

the variables A, B, C, D, E, and F will be placed in blank common in 
that order; the variables X, Y, and Z will be placed in the common area 
labeled ITEMS. 

Blank and labeled common entries appearing in COMMON statements are 
cumulative throughout the program. For example, consider the following 
two COMMON statements: 

COMMON A, B-, C /R/ D,, E /S/ F 
COMMON G, H /S/ I, J /R/P//W 

These two statements have the same effect as the single statement: 

COMMON A, B, C, G, H,, W /R/ D, E, P /S/ F, I, J 

Example 3 : 

Assume that A, B, C, K, X^, and Y each occupy four locations of 
storage, H and G each occupy eight locations, and D and E each occupy 
two locations. 

Call ing Pr ogram Subprogram 

SUBROUTINE MAPMY ( ) 

COMMON H, A /R/ X, D // B 

COMMON G, Y, C /R/ K, E 

CALL MAPMY (...) 
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Explanation ; 

In the calling program, the statement COMMON H,A/R/X,D//B causes 16 
locations (four locations each for A and B, and eight for H) to be 
reserved in blank common in the following order: 

^ ^ 

Beginning | H A B j 

of blank | 8 locations U locations 4 locations j 

common | i 

[ : ] 

I continuation of blank common I 

[ J 

and also causes six locations (four for X and two for D) to be reserved 
in the labeled common area R in the following order: 

Beginning | X D j 

of labeled! | 

common R | 4 locations 2 locations I 

, . , 

! continuation of labeled common 

I I 

L . J 

The statement COMMON G,Y,C/R/K,E appearing In ^>"^ subprogram MAPMY 

would cause the variables G, _Y — ^ <^^co share the same storctge space 

(in blank common) as H, A, tstnpT'B, respectively. It would also cause the 
variables K and id to Share the same storage space (in labeled common 
area R) as X and D, respectively. The length of a COMMON area may be 
increased by using an EQUIVALENCE statement (see the section 
"EQUIVALENCE Statements"). 

EQUI VALENCE State ment 



r -^ ' 1 

I General Form | 

|. _ ^ 

EQUIVALENCE (a„ b, c, ,..), (d, e, f , . - . ) 

Where: a, b, c, d, e, f , , . . are variables that may be subscripted. 
The suiDscripts may have two forms: If the variable is singly 
subscripted it refers to the position of the variable in the 
array (i.e., first variable, 25th variable, etc.).. If the 
variable is multi^subscripted it refers to the position in 
the array in the same fashion as the position is referred to 
in an arithmetic statement. 

The EQUIVALENCE statement provides the option for controlling the 
allocation of data storage within a single program or subprogram. It is 
analogous to the option of using the COMMON statement to control the 
allocation of data storage among several programs. In particular, when 
the logic of the program permits it, the number of storage locations 
used can be reduced by causing locations to be shared by two or more 
variables of the same or differing types and lengths. The EQUIVALENCE 
statement cannot be used to obtain mathematical equality of two 
variables. 
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Example 1 ; 

DIMENSION B(5), C(10, 10), D(5, 10, 15) 
EQUIVALENCE (A, B(l), C(5,3)), (D(5,10,2),, E) 

Explanation : 

This EQUIVALENCE Statement indicates that the variables A,B(1), and 
C(5,3) are assigned the same storage locations. In addition, it 
specifies that D(5,10,2) and E are assigned the same storage locations. 
In this case, the subscripted variables refer to the position in an 
array in the same fashion as the position is referred to in an 
arithmetic statement. Note that variables or arrays that are not 
mentioned in an EQUIVALENCE statement are assigned unique storage 
locations. The EQUIVALENCE statement must not contradict itself or any 
previously established equivalences. For example, the further equiva- 
lence specification of B{2) with any other element of the array C, other 
than C(6,3), is invalid. 

Example 2 : 

DIMENSION B(5), C(10, 10), D(5, 10, 15) 
KQUIVALENCE (A, B(l), C(25)), (D(IOO), E) 

Explanation ; 

This EQUIVALENCE stat€itn=mv~ i,«3icates that the variable A, the first 
variable in the array B, namely B(l), auci ^Ke 25th variable in the array 
C, namely C(5,3), are to be assigned the same etoraae locations. In 
addition, it also specifies that D(IOO), i.e., D(5,10,2), and E are to 
share the same storage locations. Note that the effect of the 
EQUIVALENCE Statement in examples 1 and 2 is the same. 

Variables that are brought into COMMON through EQUIVALENCE statements 
may increase the size of the block, as indicated by the following: 
statements: 

COMMON A, B, C 
DIMENSION D(3) 
EQUIVALENCE (B,D(1)) 

This would cause a common area to be established containing the 
variables A, B, and C. The EQUIVALENCE statement would then cause the 
variable D(l) to share the same storage location as B, D(2) to share the 
same storage location as C, and D(3) would extend the size of the common 
area, in the following manner: 

A (lowest location of the common area) 

B, D(l) 

C, D(2) 

D(3) (highest location of the common area) 

Since arrays must be stored in consecutive forward locations, a 
variable may not be made equivalent to another variable of an array in 
such a way as to cause the array to extend before the beginning of the 
common area. For example, the following EQUIVALENCE statement is 
invalid: 

COMMON A, B, C 
DIMENSION D(3) 
EQUIVALENCE (B, D(3)) 
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because it would force D(l) to precede A, as follows: 
D(l) 

A, D(2) Slowest location of the common area) 

B, D(3) 

C (highest location of the common area) 

Two variables in one COMMON block or in two different COMMON blocks 
may not be made equivalent. 
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SUBPROGRAMS 



It is sometimes desirable to write a program which, at various 
points, requires the same computation to be performed with different 
data for each calculation. It would simplify the writing of that 
program if the statements required to perform the desired computation 
could be written only once and then could be referred to freely, with 
each subsequent reference having the same effect as though these 
instructions were written at the point in the program where the 
reference was made. 

For example, to take the square root of a number, a program must be 
written with this object in mind. If a general program were written to 
take the square root of any number, it would be desirable to be able to 
incorporate that program (or subprogram) into other programs where 
square root calculations are required. 

The FORTRAN language provides for the above situation through the use 
of subprograms. There are three classes of subprograms: Statement 
Functions,, FUNCTION subprograms, and SUBROUTINE subprograms. In addi- 
tion, there is a group of FORTRAN supplied subprograms (see Appendix C) . 

The first two classes of subprograms are called functions. Functions 
differ from SUBROUTINE subprograms, in that functions return at least 
one value to the calling program, whereas SUBROUTINE subprograms need 
not return any. 



NAMING SUBPROGRAMS 

A subprogram name consists of from one through six alphameric 
characters, the first of which must be alphabetic (special characters 
may not be used). The type of a subprogram can be indicated in the same 
manner as variables. 

1, Typ e Dec laration of a S tatement Function ; Such declaration may be 
accomplished in one of three ways: by the predefined convention, by 
the IMPLICIT statement, or by the Explicit specification state- 
ments. Thus, the same rules for declaring the type of variables 
apply to Statement Functions. 

2- Type Dec larati on of FUNCTION Subprograms : Such declaration may be 
made in the same fashion as Statement Functions. In addition, the 
type (INTEGER, REAL, COMPLEX, and LOGICAL) may appear in the 
FUNCTION definition statement. 

3- Typ e Dec laration of a SUBROUTINE Subprogram ; The type of a 

SUBROUTINE subprogram cannot be defined, because the results that 
are returned to the calling program are dependent only on the type 
of the variable names appearing in the argument list of the calling 
program and/or the implicit arguments in COMMON. 



FUNCTIONS 

A function is a statement of the relationship between a number of 
variables. To use a function in FORTRAN, it is necessary to: 
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Define the 
performed) , 



function (i.e., specify what calculations are to be 



2, Refer to the function by name, where required in the program. 



Function Definitio n 

There are three steps in the definition of a function in FORTRAN: 

1, The function must be assigned a unique name by which it may be 
called (see the section "Naming Subprograms") . 

2, The arguments of the function must be stated. 

3, The procedure for evaluating the function must be stated. 

Items 2 and 3 are discussed in detail in the sections dealing with 
the specific subprogram (e.g., "Statement Functions," "FUNCTION Subpro- 
grams , " etc . ) . 



Function Reference 



When the name of a function appears in any FORTRAN arithmetic 
expression, this, effectively, references the function. Thus,, the 
appearance of a function with its arguments in parentheses causes the 
computations to be performed as indicated by the function definition- 
The resulting quantity replaces the function reference in the expression 
and assumes the type of the function. The type and length of the name 
used for the reference must agree with the type and length of the name 
used in the definition. 



STATEMENT FUNCTIONS 



Statement functions are defined by a single arithmetic or logical 
assignment statement within the program in which they appear. For 
example, the statement: 

FUNC(A,B) = 3,*A+B**2.+X+Y+Z 

defines the statement function FUNC, where FUNC is the function name and 
A and B are the function arguments. 

The expression on the right defines those computations which are to 
be performed when the function is used in an arithmetic statement. This 
function might be used in a statement as follows: 

C = FUNC(D,E) 

which is equivalent to: 

C = 3,*D+E**2.+X+Y+Z 

Note the correspondence between A and B in the function definition 
statement and D and E in the arithmetic statement. The quantities A and 
B enclosed in parentheses following the function name are the arguments 
of the function. They are dummy variables for which the quantities D 
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and E, respectively are substituted when the function is used in an 
arithmetic statement. 



r 1 

I General Form | 

^ „_^ 

name (a,b, ...»n) = expression 

Where: name is any subprogram name (see the section "Naming 
Subprograms" ) , 

a,b, --.»n are distinct (within the same statement) nonsub- 
scripted variables, 

expression is any arithmetic or logical expression that does 
not contain subscripted variables. Any statement functions 
appearing in this expression must be defined previously. 

L , ^ J 

A maximum of 15 variables appearing in the expression may be used as 
arguments of the function. 

Note ; All Statement Function definitions to be used in a program must 
precede the first executable statement of the program. 

Examples ; 

Valid statement function definitions: 



SUM(A,B,C,D) = A+B+C+D 
FUNC(Z) = A+X*Y*Z 
AVG(A,B,C,,D) = (A+B+C+D) /4 
ROOT(A,B,C) = SQRT(A**2+B**2+C**2) 
VALID (A, B) = , NOT. A, OR, B 

Note: The same dummy arguments may be used in more than one Statement 
Function definition and as variables outside Statement Function defini- 
tions. 

Invalid statement function definitions: 



SUBPRGO, J,K)=3*I+J**3 (arguments must be variables) 
SOMEF(A(I) ,B)=A(I)/B+3, (arguments must be nonsub- 

scripted) 
SUBPROGRAM(A, B)=A**2+B**2 (function nam.e exceeds limit 

of six characters) 
3FUNC(D)=3.1U*E (function name must begin with 

an alphabetic character) 
ASF(A)=A+B(I) (subscripted variable in the 

expression) 

Valid statement function references : 

NET = GROS - SUM(TAX,, FICA, HOSP,, MISC) 

ANS = FUNC( RESULT) 

GRADE = AVG(LAB, LECTUR, SUM(TEST1, TEST2, TEST3 , TEST4), FACTOID) 
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Invalid statement function references ; 

WRONG = SUM (T.?^, PICA) (number of arguments 

does not agree with 
above definition) 

MIX = FUNCCI) (mode of argument 

does not agree with 
above definition) 



FUNCTION SUBPROGR.AMS 

The FUNCTION subprogram is a FORTRAN subprogram consisting of any 
number of statements. It is an independently written program that is 
executed wherever its name appears in another program. 

r 1 

1 General Form | 

^ _ ^ 

FUNC TI ON name ( a j. , a 2 « a 3 » - • • » an ) 



RETURN 



END 

Where: name is subprogram name (see the section "Naming 
Subprograms" ) , 

a.i.»a2#a3„ . . . ,an are nonsubscripted variable, array, or dummy 
names of SUBROUTINE or other FUNCTION subprograms. (There 
must be at least one argument in the argument list. ) 

L . J 

Since the FUNCTION is a separate subprogram, the variables and 
statement numbers within it do not relate to any other program- 

The FUNCTION subprogram may contain any FORTRAN statement except a 
SUBROUTINE statement, another FUNCTION statement,, or BLOCK DATA state- 
ment. 

The arguments of the FUNCTION subprogram (i.e., §.1,32, as , . . . ,an) may 
be considered to be dummy variable names. These are replaced at the 
time of execution by the actual arguments supplied in the function 
reference in the calling program. The actual arguments may be any of 
the following: any type of constant,, any type of subscripted or 
nonsubscripted variable, an arithmetic or logical expression, or the 
name of another subprogram- The actual arguments must correspond in 
number, order, and type to the dummy arguments. The array size must 
also be the sarae» except when adjustable dimensions are used. All 
arguments in a subprogram refer to the storage area assigned to the 
arguments by the calling program. 

The relationship between variable names used as arguments in the 
calling program and the dummy variables used as arguments in the 
FUNCTION subprogram is illustrated in the following example: 
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Example 1 : 

Calling Program 

A = SOMEF(B,C) 



FUNCTION Subprogram 

FUNCTION SOMEF(X,Y) 
SOMEF = X/Y 
RETURN 
END 



Explanation ; 

In the above example, the value of the variable B of the calling 
program is used in the subprogram as the value of the dummy variable X; 
the value of C is used in place of the dummy variable Y. Thus if B = 
10.0 and C = 5.0, then A = B/C, which is equal to 2.0. 

The name of the function must be assigned a value at least once in 

the subprogram as the argument of a CALL statement, as a DO variable, as 

the variable name on the left side of an arithmetic statement, or in an 
input list (READ statement) within the subprogram. 



Example 2 ; 

Calling Program 



FUNCTION Subprogram 
FUNCTION CALC (A,B,J) 



I = J*2 
ANS = R00T1*CALC(X,Y,I) 



CALC = A**I/B 



RETURN 
END 

Explanation : 

In this example, the values of X, Y, and I are used in the FUNCTION 
subprogram as the values of A, B, and J, respectively. The value of 
CALC is computed, and this value is returned to the calling program, 
where the value of ANS is computed. The variable I in the argument list 
of CALC in the calling program is not the same as the variable I 
appearing in the subprogram. 

When a dummy argument is an array name, an appropriate DIMENSION or 
Explicit specification statement must appear in the FUNCTION subprogram. 
None of the dummy arguments may appear in an EQUIVALENCE statement or a 
COMMON statement. 



Type Specification of the FUNCTI ON Subprogram 



In addition to the three ways of declaring the type of a FUNCTION 
name (i.e., predefined conveintion, IMPLICIT statement. Explicit specifi- 
cation statement) , there exists the option explicitly specifying the 
type of a FUNCTION name within the FUNCTION statement. 
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r ■ ' T 

I General Form | 

^ ^ 

Type FUNCTION name*s (aj., a2,a3, . . - ,,an) 

Where: Type is INTEGER, REAL, COMPLEX, or LOGICAL. 

name is the name of the FUNCTION subprogram. 

*s is optional and represents one of the permissible length 
specifications for its associated type. 

cia.#a2,a3, .«- ,an are nonsubscripted variable,, array, or dummy 
names of SUBROUTINE or other FUNCTION subprograms. (There 
must be at lea^t one argument in the argument list. ) 

L . J 



Example 1 ; 

REAL FUNCTION SOMEF (A, B) 

SOMEF = A**2 + B**2 

RETURN 
END 

Example 2 : 

INTEGER FUNCTION CALC*2 (X,Y,,Z) 

CALC = X+Y+Z**2 



RETURN 
END 

Explanation : 

The FUNCTION subprograms SOMEF and CALC in examples 1 and 2 are 
declared as type REAL (length 4) and INTEGER (length 2) , respectively. 



RETURN and END Statements in a Function Subprogram 

Note that all of the preceding examples of FUNCTION subprograms 
contain both an END and at least one RETURN statement. The END 
statement specifies, for the compiler, the end of the subprogram; the 
RETURN statement signifies a logical conclusion of the computation and 
returns any computed value and control to the calling program. 

There may, in fact, be more than one RETURN statement in a FORTRAN 
subprogram, 

S ubprog rams 8 5 



Example ; 



FUNCTION DAV (D,E,F) 
IF (D-E) 10, 20^ 30 
10 A = D+2.0+E 



5 A = F+2,0*E 



20 DAV = A+B**2 



RETURN 
30 DAV = B**2 



RETURN 
END 



SUBROUTINE SUBPROGRAMS 



The SUBROUTINE subprogram is similar to the FUNCTION subprogram in 
many respects: the rules for naming FUNCTION and SUBROUTINE subprograms 
are the same, they both require an END statement, and they both contain 
the same sort of dummy arguments. Like the FUNCTION subprogram, the 
SUBROUTINE subprogram is a set of commonly used computations, but it 
need not return any results to the calling program, as does the FUNCTION 
subprogram. 

The SUBROUTINE subprogram is called by the CALL statement, which 
consists of the word CALL followed by the name of the subprogram and its 
parenthesized arguments. 



j General Form 
\- 



SUBROUTINE name (3.1,3.2,33 ,,.. ,an) 

RETURN 

END 

where: name is the subprogram name (see the section "Naming 
Subprograms" ) . 

^±»sL2r^3* ' • • i^n are arguments. (There need not be 
any.) ""Each argument used must be a nonsubscripted variable 
or array name, the dummy name of another SUBROUTINE or 
FUNCTION subprogram, or of the form * where the character 
"*" denotes a return point specified by a statement number 
in the calling program. 



Since the SUBROUTINE is a separate subprogram, the variables and 
statement numbers within it do not relate to any other program. 
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The SUBROUTINE subprogram may contain any FORTRAN statememt except a 
FUNCTION statement, another SUBROUTINE statement, or BLOCK DATA state- 
ment - 

The SUBROUTINE subprogram may use one or more of its cirguments to 
return values to the calling program. Any arguments so used must appear 
on, the left side of an arithmetic statement or in an input list within 
the subprogram, as arguments of a CALL statement, or as DO variables. 
The SUBROUTINE name must not appear in any other statement in the 
SUBROUTINE subprogram. 

The arguments (a^., a.^^ a3,-..,an) may be considered duirany variable 
names that are replaced at the time of execution by the actual arguments 
supplied in the CALL statement. The actual arguments must correspond in 
number, order, and type to the dummy arguments. The array size must 
also be the same except when adjustable dimensions are used. Dummy 
arguments may not appear in an EQUIVALENCE or DATA statement within the 
subprogram nor may they be given initial data values in an Explicit 
specification statement. 

Example ; The relationship between variable names used as arguments in 
the calling program and the dummy variable used as arguments in the 
SUBROUTINE subprogram is illustrated in the following example. The 
object of the subprogram is to "copy" one array directly into another- 

Main Program SUBROUTINE Subprogram 

DIMENSION X(100),Y(100) 

SUBROUTINE COPY(A,B,N) 

DIMENSION A (100),B(100) 

DO 10 I = 1, N 
CALL COPY (X, Y, K) 10 B(I) = A (I) 

RETURN 

END 



CALL Statement 

The CALL statement is used only to call a subroutine subprogram- 

r ' 1 

I General Form | 

|.__ _ ^ 

CALL name ( a^ , aa , 3.3 , . . . , an) 

Where: name is the name of a subroutine subprogram. 

§.i.»S.2»§.3# •- • #£n ^^^ the actual arguments that are being 
supplied to the subroutine subprogram. Each may be of the 
form Bn where n is a statement number (see "RETURN State- 
ments in a SUBROUTINE Subprogram" ) . 

L J 

Examples ; 

CALL OUT 

CALL MATMPY (X, 5, 40, Y, 7, 2) 
CALL QDRTIC (X, Y, Z, R00Tl„R00T2) 
CALL SUBl {X+Y*5, 'ABDF* ,SINE) 
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The CALL statement transfers control to the subroutine subprogram and 
replaces the dummy variables with the value of the actual arguments that 
appear in the CALL statement. The arguments in a CALL statement may be 
any of the following: any type of constant, any type of subscripted or 
nonsubscripted variable, an arithmetic expression,, the name of a 
subprogram, or a statement number (see "RETURN Statements in a 
SUBROUTINE Subprogram" ) . 

The arguments in a CALL statement must agree in number, order,, and 
type with the corresponding arguments in the subroutine subprogram. The 
array sizes must also be the same in the subroutine and the calling 
programs, except when adjustable dimensions are used (see "Adjustable 
Dimensions"). If an actual argument corresponds to a dummy argument 
that is defined or redefined in the referenced subprogram, the actual 
argument must be a variable name, subscripted variable name, or array 
name. All arguments in a subprogram refer to the storage area assigned 
to the arguments by the calling program. 



RETURN Statement in a SUBROUTINE Subprogram 



r 1 

I General Form | 

|._^ -^ 

RETURN 

RETURN i 

Where: i is an integer constant or variable of length 4 whose value, 
say n, denotes the nth statement number in the argument list 
of a SUBROUTINE statement. 



The normal sequence of execution following the RETURN statement of a 
SUBROUTINE subprogram is to the next statement following the CALL in the 
calling program. It is also possible to return to any numbered 
statement in the calling program by using a return of the type where i 
is an integer constant or variable. Returns of the type RETURN may be 
made in either a SUBROUTINE or FUNCTION subprogram (see, "RETURN and END 
Statements in a FUNCTION Subprogram"). Returns of the type RETURN i may 
only be made in a SUBROUTINE subprogram. In a main program, a RETURN 
statement performs the same function as a STOP statement. 
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Example ; 

Calling Progra m Subprogram 



SUBROUTINE SUB (X, Y, Z, *,, ♦) 



10 CALL SUB (A,B,C,, &30, gUO) 

20 Y = A + B 100 IF (R) 200,300,400 

200 RETURN 

300 RETURN 1 

400 RETURN 2 

30 Y = A + C END 



40 Y = B + C 



END 

Explanation ; 

In the preceding example, execution of statement 10 in the calling 
program causes entry into subprogram SUB, When stateiment 100 is 
executed, the return to the calling program will be to stateiment 20, 30, 
or 40,, if R is less than, equal to, or greater than zero, respectively. 

A CALL statement that uses a RETURN i form may be best understood by 
comparing it to a CALL and computed GO TO statement in sequence. For 
example, the CALL statement; 

CALL SUB (P,&20,,Q,€35,R,S22) 

is equivalent to; 

CALL SUB (P,Q,,R.I) 
GO TO (20,35, 22), I 

where the index I is assigned a value of 1, 2, or 3 in the called 
subprogram- 



Multiple ENTRY into a Subprogram 

The standard (normal) entry into a SUBROUTINE subprogram from the 
calling program is made by a CALL statement that references the 
subprogram name. The standard entry into a FUNCTION subprogram is made 
by a function reference in an arithmetic expression. Entry is made at 
the first executable statement following the SUBROUTINE or FUNCTION 
statement. 

It is also possible to enter a subprogram (either SUBROUTINE or 
FUNCTION) by a CALL statement or a function reference that references an 
ENTRY statement in the subprogram. Entry is made at the first 
executable statement following the ENTRY statement. 
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j General Form 

^ 

ENTRY name (a^f §,2,33, . . . ,an) 

Where: name is the name of an entry point (see "Naming 
Subprograms" ) . 

3.^,32,33, ., .,an are the dummy arguments corresponding to an 
actual argument in a CALL statement or in a function 
reference. 



ENTRY statements do not affect control sequencing during normal 
execution of a subprogram- The order, type, and number of arguments 
need not agree between the SUBROUTINE or FUNCTION statement and the 
ENTRY statements, nor do the ENTRY statements have to agree among 
themselves in these respects. Each CALL or function reference,, however, 
must agree in order, type, and number with the SUBROUTINE, FUNCTION, or 
ENTRY statement that it references. Entry may not be made into the 
range of a DO; further, a subprogram may not reference itself directly 
or through any of its entry points. This statement is regarded as 
nonexecutable within its subprogram. If it appears in a function 
subprogram the name given in the FUNCTION statement is still used to 
return the value of the function to the point of reference, rather than 
the name of the ENTRY statement. 



Example 1 ; 

Calling Program 



Subprogram 

SUBROUTINE SUBl (U, V, W, X, Y, Z) 



1 CALL SUBl (A, B,C,D»E,F) 



U = V 



2 CALL SUB2 (G, H,P) 



ENTRY SUB 2 (T,U,V) 



3 CALL SUB 3 



ENTRY SUB 3 



END 



Explanation ; 

In the preceding example, the execution of statement 1 causes entry 
into SUBl, starting with the first executable statement of the 
subroutine- Execution of statements 2 and 3 also causes entry into the 
called program, starting with the first executable statement following 
the ENTRY SUB2(T,U,V) and ENTRY SUB3 statements, respectively. 

Entry into a subprogram initializes all references in the whole 
subprogram to items in the argument list. Return from a subprogram is 
made by way of the entry point referenced. ENTRY statements may only 
appear in FUNCTION or SUBROUTINE subprograms. The dummy arguments in a 
subprogram may appear in any statement if they first appear as dummy 
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arguments in a FUNCTION, SUBROUTINE, or ENTRY statement. The following 
is a valid example: 

SUBROUTINE SUB (X,Y,Z,I) 



ENTRY SUBl (A, B) 

C = A+B 

Example 2 ; 

Calling Program Subprogram 

SUBROUTINE SUBl (U, V, W, X, Y, Z) 

RETURN 

ENTRY SUB 2 (T, *,♦) 

U = V* W+T 

ENTRY SUB3 (♦,♦) 

X = Y**Z 

IF (U-X) 100, 200, 300 

RETURN 1 

RETURN 2 

RETURN 

END 



CALL SUBl (A,B,C,D,E,F) 



CALL SUB2(G,S10, g20) 



50 
100 
200 
300 



CALL SUB3(g 10^8 20) 
5 Y =A+B 



10 Y = C+D 
20 Y = E+F 



Explanation ; 

In the example above, a call to SUBl merely performs initialization. 
Subsequent calls to SUB2 and SUB3 result in execution of different 
sections of the subroutine SUBl. Then, depending upon the result of the 
arithmetic IF at statement 50, return is made to the calling program at 
statement 10, 20, or 5. 



Additional Rules for Using ENTRY 



A CALL may only change the value of explicit arguments (or implicit 
arguments in COMMON). It cannot affect the value of those which 
are initialized by some previous CALL. 

If a name is identified as a dummy argument only by its appearance 
in a given ENTRY statement, no use of that dummy argument may 
appear in statements preceding (physically) the ENTRY statement. 
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3. The appearance of an ENTRY statement does not alter the rules 
regarding the placement of Statement Functions in subprograms. 



The EXTERNAL Statement 



I General Form j 

|. _ ^ 

EXTERNAL a^b,c, 

Where: a,b,,c, .,, are names of subprograms that are used as 
arguments in other subprograms. 

L J 



The name of any subprogram that is used as an argument in another 
subprogram must appear in an EXTERNAL statement. For example, assume 
that SUB and MULT are subprogram names in the following statements : 



Example 1 : 

Calling Program 

EXTERNAL MULT 



Subprogram 

SUBROUTINE SUB (X, Y, Z) 
IF (X) 4,6,6 
D = Y (X,Z**2) 



CALL SUB (A, MULT,C) 



RETURN 
END 



Explanation ; 

In this example, the subprogram name MULT is used as an argument in 
the subprogram SUB,. The subprogram name MULT is passed to the dummy 
variable Y, as are the variables A and C passed to the dummy variables X 
and Z, respectively. The subprogram MULT will be called and executed 
only if the value of A is negative. 

Example 2 



CALL SUB (A, B, MULT (C,D),37) 



Explanation ; 

In this example, an EXTERNAL statement is not required because the 
subprogram named MULT is not an argument; it is executed first and the 
result becomes the argument. 
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FORTRAN SUPPLIED SUBPROGRAMS 

FORTRAN provides the programmer with a library of commonly used 
function and subroutine subprograms. There are three classes of 
subprograms provided: 

1, Mathematical function subprograms — defined as FUNCTION subprograms. 

2^ Subroutines which test the status of pseudo machine indicators 
(sense lights) — defined as SUBROUTINE subprograms. 

3. The three subprograms EXIT,, DUMP, and PDUMP — also defined as 
SUBROUTINE subprograms. 

The EXIT subroutine terminates program execution; DUMP dumps storage 
and terminates program execution; PDUMP dumps storage and continues 
program execution. 

Variables used as arguments of mathematical functions must be 
declared (i.e.* by the Explicit specification statements, the IMPLICIT 
statement, or the predefined convention) , in accordance with the 
function in which they appear. 

The entire library^ along with appropriate definitions of each 
subprogram, is given in Appendix C. 



BLOCK DATA SUBPROGRAM 

In order to enter data into a COMMON block, a separate subprogram 
must be written. This separate subprogram contains only the DATA, 
COMMON, DIMENSION, EQUIVALENCE, and Type statements associated with the 
data being defined. Data may be entered into labeled (named), but not 
unlabeled, COMMON by the BLOCK DATA subprogram. 

r ' ' ^ 1 

I General Form | 

^ „ ^ 

BLOCK DATA 



END 

L . : J 

1, The BLOCK DATA subprogram may not contain any executable state- 
ments. 

2. The first statement of this subprogram must be the BLOCK DATA 
statement. 

3,. All elements of a COMMON block must be listed in the COMMON 
statement, even though they do not all appear in the DATA 
statement; for example, the variable A in the COMMON statement in 
the following example does not appear in the DATA statement: 

BLOCK DATA 

COMMON/ELN/C, A, B/RMG/Z, Y 

REAL B(U)/1.0,1. 2,2*1.3/,Z*8(3)/3*7.6U980825D0/ 

COMPLEX C/ (2. 4*3, 769)/ 

END 
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4, Data may be entered into more than one COMMON block in a single 
BLOCK DATA subprogram. 

5- No element may have more than one initial value assigned in the 
same program- 
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APPENDIX A: SOURCE PR0GRAI4 CHARACTERS 



r- 






-T- 


■"— ■ 


""^ 






-T" 


»_«..._«».^ . 


-T- 


. — . _— — _^ 


■" 


1 




Alphabetic 




EB( 


2DIC or 


BCDIC 




Numeric 




EBCDIC 


or 


BCDIC 1 




Characters 




Card 


Punches 




Characters 




Card Punches | 


L_ 






1 










-4-- 




-4— 






J 


r 


A 
B 
C 
D 
E 
F 
G 
H 
I 
J 
K 
L 
M 




t 






12-1 
12-2 
12-3 
12-4 
12-5 
12-6 
12-7 
12-8 
12-9 
11-1 
11-2 
11-3 
11-4 




t 
1" 



1 
2 
3 
4 
5 
6 
7 
8 
9 


T 
-J— 





1 
2 
3 
4 
5 
6 
7 
8 
9 


1 

J 






T 




— i- 


1 




N 










11-5 






Special 




EBCDIC 


1 


BCDIC 1 















11-6 






Characters 


1 Card Punches | Card Punches | 




P 










11-7 




1-- 




-4— 




_1 


4 






T 




T 


1 




Q 










11-8 






+ 




12-6-8 




12 1 




R 










11-9 






- 




11 




11 1 




S 










0-2 






/ 




0-1 




0-1 1 




T 










0-3 






= 




6-8 




3-8 1 




U 










0-4 






1, 




12-3-8 




12-3-8 1 




V 










0-5 






) 




11-5-8 




12-4-8 1 




w 










0-6 






* 




11-4-8 




11-4-8 1 




X 










0-7 






, (comma) 


1 


0-3-8 




0-3-8 1 




Y 










0-8 






( 




12-5-8 




0-4-8 1 




Z 










0-9 




1 


(apostrophe) 


1 


5-8 




4-8 1 




$ 










11-3- 


8 




blank 


1 


[no punch) 


1 (no punch) | 


I. 






_J_^ 










_X_ 




_j._ 




_x 


J 


r 


Source 


progra 


ms 


are coded in 


either BCDIC or 


EBCDIC charact 


1 

-er codes. | 




Mixing 


of 1 


the 


two 


, however,. 


is 


not allowed. 








1 


L. 


























J 



The 4 8 characters listed above comprise the set of characters 
acceptable by FORTRAN. In previously implemented FORTRAN languages, 
there existed dual characters in the sense that two graphics 
(characters) were associated with a single card code. The most commonly 
used set of dual characters was the following: 



(apos, ) 



However, in IBM System/360 each of these duals now has sejparate card 
codes. Thus, when specifying, for instance, a + character^ care should 
be taken that a 12-6-8 punch is used instead of a 12 punch which now 
represents an 6. The card codes for the remaining duals is as follows: 
#(8-3), (this graphic is now represented as a 0(12-4-8), % (0-4-8), 
and a (8-4). 
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APPENDIX B; OTHER FORTRAN FEATURES ACCEPTED BY FORTRAN IV 



This section discusses those features of previously implemented 
FORTRAN IV languages that are incorporated into the IBM Time Sharing 
System/360 FORTRAN IV language. The inclusion of these additional 
language facilities allows existing FORTRAN programs to be re-compiled 
for use on the IBM System/360 with little or no re- prog ramming. 



READ Statement 



j General Form 
h 



READ b, list 



Where: 



b, is the statement number or array 
statement describing the data. 



name of the FORMAT 



list is a series of variable or array names, separated by 
commas, which may be indexed and incremented- They specify 
the number of items to be read and the locations in storage 



into which the data is placed. 



This statement causes data to be read from the 
with the system input. 



data set associated 



PUNCH Statement 



j General Form 
k 



PUNCH b, list 



Where: 



b is the statement number or 
statement describing the data. 



array name of the FORMAT 



list is a series of variable or array names, 
commas^ which may be indexed and incremented, 
the number of items to be written and the 
storage from which the data is taken. 



separated by 

They specify 

locations in 



The PUNCH statement causes data to be written in the data set 
associated with the system output. 
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PRINT Statement 



J. . . ^ . ^ 

I General Form | 

|. _ _ ^ ^ 

PRINT b. list 

Where: b is the statement number or array name of the FORMAT 
statement describing the data, 

list is a series of variable or array names, separated by 
commas, which may be indexed and incremented. They specify 
the number of items to be written and the locations in 
storage from which the data is taken. 

L . , , J 

The PRINT statement causes data to be written in the; data set 
associated with the system output. 



DATA Initialization Statement 



r ■ 1 

I General Form | 

^ ^ ^ 

DATA Vi,..,.,Vn/ii*di,...,in*dn/»Vn+i.,...»v /in+i*dn+i, - • • ,i *d /,... 

Where: yi,,.,,,v are variables,, subscripted variables (in which 
case the subscripts must be integer constants),, or array 
names . 

di,.,,,d are values representing integer, real, complex, 
logical, or literal data constants. 

ii,-.-»i represent unsigned integer constants indicating 
the nuniber of consecutive variables that are to be assigned 
the value of di,,..,d . 

L Z. u . J 

A data initialization statement is used to define initial values of 
variables and arrays. There must be a one-for-one correspondence 
between these variables (i.e., Vi,*..,v ) and the data constants (i.e., 
di. , . . . y d ) • 

Example 1 : 

DIMENSION D(5,10) 

DATA A, B, C/5. 0,6. 1,7, 3/, 0/25*1. 0/ 

Explanation ; 

The DATA statement indicates that the variables A„ B, and C are to be 
initialized to the values 5.0, 6,1, and 7,3, respectively. In addition, 
the statement specifies that the first 25 variables in the array D are 
to be initialized to the value 1.0. 

Example 2 : 

DIMENSION A(5), B(3,3), L(U) 

DATA A/5*l,0/^ B/9*2,0/, L/U + .TRUE./, C/'FOURV 
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Explanation ; 

The DATA statement specifies that all the variables in the arrays A 
and B are to be initialized to the values 1.0 and 2.0, respectively. 
All the logical variables in the array L are initialized to the value 
.TRUE-. The letters T and F may be used as an abbreviation for .TRUE- 
and .FALSE.,, respectively. In addition, the variable C is initialized 
with the literal data constant FOUR. 

An initially defined variable, or variable of an array, may not be in 
blank common; however, in a labeled common block, they may be initially 
defined only in a block data subprogram. (See the section 
"SUBPROGRAMS.") 



DOUBLE PRECISION Statement 



r 1 

I General Form | 

^ __^ 

DOUBLE PRECISION a,b,c, . . . 

Where: a,b#c,..- are variable names that may be dimensioned in the 
statement, or function names. 

L ^ —J 

The DOUBLE PRECISION statement explicitly specifies that the vari- 
ables a,b,c, ... are of type double precision. This statement overrides 
any specification of a variable made by either the predefined convention 
or the IMPLICIT statement. This specification is identical to that of 
type REAL* 8. 

In addition, FUNCTION subprograms may be typed double precision,, as 
follows: 

DOUBLE PRECISION FUNCTION name (a^., 32,33, ... ,an) 



Arguments of a FUNCTION or SUBROUTINE Program Enclosed by Slashes 

Arguments in a FUNCTION or SUBROUTINE subprogram may be enclosed in 
slashes within the commas. This form is equivalent to the normal format 
without the slashes. 
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APPENDIX C: 



FORTRAN SUPPLIED SUBPROGRAMS 



MATHEMATICAL FUNCTION SUBPROGRAMS 



Table 4. Mathematical Function Subprograms 

r T T 
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1 Definition 
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j Comp] 
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1 CDEXP 
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JReal 
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|DLOG 
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|Real 
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*8 
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.ex 


*8 




1 CDLOG 
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j Comp] 
_X 


ex 


+ 16 
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_x 
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(Continued) 
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Table 4. Mathematical Function Subprograms (Continued) 
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JReal *4 
JReal *8 
-X 


+u 


1 Integer 
JReal *4 
JReal *8 
X 


♦ 4 






T - 


T 


— r 




T 




T 




T 






1 Absolute value 


|IABS 

1 ABS 

|DABS 
_ 1 


llArgI 


-X— 


I 
I 
I 


X- 


1 
1 

1 


1 Integer 
IReal +4 
IReal +8 

-X 


*4 


1 Integer 
JReal +4 
IReal *8 
X 


♦ 4 






-f 




T- 




1 




t 




t 






1 Modulus 


|CABS 


1 J a2+b2 1 for a+bi 









1 


1 Complex 


♦ 8 


IReal +4 








ICDABS 


X — 


_4. 





-X 


1 


1 Complex 
. J. 


♦ 16 


|Real +8 

-X 








T — 


t 


t 




1 




t — 




t 






[Truncation 


|INT 
|AINT 


1 Sign of Arg times 
[largest integer 
|<|Arg| 




I 
I 




1 
1 


IReal +4 
IReal *4 




1 Integer 
IReal *4 


♦ 4 






1 IDINT 
-4- 


-4- - - 


_x 


I 


1 


1 


IReal ♦S 

.X 




1 Integer 
X - - 


*4 






T 


t 


T 




T 




T 




T 






1 Largest value 


lAMAXO 
JAMAXl 


jMax (Argi, Argar • . . 


) 1 


I 
I 




>2 
>2 


1 Integer 
JReal +4 


♦ 4 


IReal *4 
IReal *4 








|MAXO 






I 




>2 


1 Integer 


♦ 4 


1 Integer 


♦ 4 






|MAX1 






I 




>2 


IReal *4 




1 Integer 


+ 4 






IDMAXI 
L 


_ 1 


X 


I 


X- 


>2 


IReal +8 
_j. 




IReal *8 
X 








T — 


T 


T 




T 




t 




T 






[Smallest value 
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JAMINl 


|Min (Argi, Arga, . - . 


) 1 


I 
I 




>2 
>2 
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JReal *4 


♦ 4 


IReal *4 
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|MINO 






I 




>2 


1 Integer 


*4 


1 Integer 


♦ 4 






IMINl 






I 




>2 


IReal +4 




1 Integer 


♦ 4 






IDMINI 


-J- — 


_x 


I 


1 


>2 


IReal *8 
-X 




|Real *8 
„x 








t 


T 


T 




1 




T 




T 






1 Float 


1 FLOAT 


1 Convert from 




I 




1 


1 Integer 


+ 4 


IReal *4 








|DFLOAT| integer to real 




I 




1 


1 Integer 


*4 


IReal *8 








1 


X _ _ 


1 




J 




1 




X 








t 


-+ - - 


t 




1 




1 




T 






jFix 


|IFIX 


1 Convert from 




I 




1 


IReal *4 




1 Integer 


♦ 4 






|HFIX 
_4. 


|real to integer 
_x _ _ 


-X 


I 


I 


1 


IReal *4 
-X 




1 Integer 
X 


*2 






t 


T - - 


T 




T 




T 




T 






1 Transfer of sign 


[SIGN 


1 Sign of Arga times 
1 1 Argi 1 




I 




2 


IReal *4 




IReal *4 








IISIGN 






I 




2 


1 Integer 


*4 


1 Integer 


*4 






IDSIGN 
1 _ 


I 


1 


I 


!_ 


2 


IReal *8 

-X 




IReal +8 








-t 


t 


T 




1 




T 




1 






1 Positive 


|DIM 


1 Arg^-Min ( Arg^. , Arg^ 


) 1 


I 




2 


IReal +4 




IReal +4 






1 difference 


|IDIM 
4. 


_4. 


X 




1 




1 Integer 
-X 


♦ 4 


1 Integer 
X 


*4 






T 


t 


T 




T 




T 




T 






1 Obtaining most 


|SNGL 






I 




1 


IReal *8 




IReal *4 






1 significant part 
























|of a Real *8 
























1 argument 


_X 


-X — 


X 




1 




-X 




-X 








T 


t 


T 




T 




T 




T 






(Obtain real 


IREAL 






I 




1 


1 Complex 


+ 8 


IReal +4 






1 part of complex 
























1 argument 
























L _ 


_i 


_X 


_X_ 




_X_ 




_x 





_x 





_J 



(Continued) 
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Table U. Mathematical Function Subprograms (Continued) 



Function 



Name 



Definition 



i In-Line (I) 
jout-of-Line (O) 



No. of 
Arg. 



Type of 
Arguments 



■f- 



Function 



Obtain imaginary 
part of complex 
argument 



AIMAG 



Complex +8 



Real +4 



Express a Real 
*4 argument in 
Real *8 form 



DBLE 



Real *4 



Real *8 



Express two real 
arguments in com- 
plex form 



CMPLX 
DCMPLX 



C=Argi,+iArg2 



Real *4 
Real *8 



Complex *8 
Complex *16 



Obtain conjugate 
of a complex 
argument 



CONJG 
DCONJG 



C=X-iY 

For Arg=X+iY 



Complex +8 
Complex *16 



Complex *8 
Complex *16 
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MACHINE INDICATOR TESTS 

In the following list of pseudo machine indicator test subroutines, 
assume that i is an integer expression and that j is an integer 
variable. These subroutines are referred to by CALL statements. 

SLITE (i) ; If i = 0, all sense lights will be turned off. If i = 1, 2, 
3,, or 4, the corresponding sense light will be turned on. 



SLITET (i> j) : Sense light i (equal to 1, 2, 3, or 4) will be tested and 
turned off. The variable j will be set to 1 if i was on, or j will be 
set to 2 if i was off. 

Example : 

Assume that it is desired to continue with the program if sense light 
i is on and to write results if sense light i is off. This can be done 
by using the Logical IF statement or a computed GO TO statement: 



CALL SLITET (3, KEN) 

GO TO (6, 17) ,KEN 
17 WRITE (3, 26) (ANS (K) , K=l, 10) 
6 CONTINUE 



Explanation: 

When the statement CALL SLITET (3, KEN) is executed, the variable KEN 
is assigned the value 1 or 2 depending on whether sense light 3 is on or 
off, respectively (and the sense light is turned off) . If KEN is 1,, 
statement 6 is executed next; if KEN is 2, statement 17 is executed. 

OVERFL ( i ) : j is set to 1 if a floating-point overflow condition 
exists, i.e., if the result of an arithmetic operation is greater than 
16^3, • j is set to 2 if no overflow condition exists; j is set to 3 if 
floating-point underflow condition exists, i.e., if the result of an 
arithmetic operation is less than 16-®^. The machine is left in a no 
overflow condition. 

DVCHK ( i ) : If the divide check indicator is on, j is set to 1 and the 
divide check indicator is turned off; if the divide check indicator is 
off, j is set to 2. 



THE EXIT> DUMP, AND PDUMP SUBPROGRAMS 

EXIT Subprogram 

A CALL to the EXIT subprogram terminates the execution of the object 
program. 
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DUMP Subprogram 

A CALL to the DUMP subprogram by the statement: 

CALL DUMP (A3L,Bi,Fi,..,,,An,rBn»Fn) 

causes the indicated limits of storage to be dumped and execution to be 
terminated, 

1. A and B are variable data names that indicate the limits of storage 
to be dumped? either A or B may represent upper or loweir limits. 

2. Fn is an integer indicating the dump format desired: 



Fn = 


Hexadecimal 


1 


Logical *1 


2 


Logical *4 


3 


Integer *2 


4 


Integer *H 


5 


Real *i\ 


6 


Real *8 


7 


Complex *8 


8 


Complex* 16 


9 


Literal 



If the argument Fn is omitted, it is assumed to be equal to 0,, and 
the dump will be hexadecimal. 

The arguments A and B should be in the same program (main program 
or subprogram) or same common block. 



PDUMP Subprogram 

A CALL to the PDUMP subprogram by the statement: 

CALL PDUMP (Ai, Bj., F^.,, . . . , An»Bn,, Fn) 

causes the indicated limits of storage to be dumped and execution to be 
continued. The PDUMP arguments are the same as the DUMP arguments. 
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APPENDIX D: SAMPLE PROGRAMS 



SAMPLE PROGRAM 1 



Sample program 1 (Figure 2) is designed to find all of the prime 
numbers between 1 and 1000. A prime number is an integer that cannot be 
evenly divided by any integer except itself and 1. Thus 1, 2, 3, 5, 7, 
11,.,. are prime numbers. The number 9, for example, is not a prime 
number, since it can evenly be divided by 3. 



IBM FORTRAN Codlnii Porm "'.'""l 


'"— " SAMPLE PROGRAM 1 


:T,i^^,=o.s 


G..,«,C 
















„« 1 o, 1 


— 1°- 6/66 


PUNCH 


















"»™«' 


^ FORTRAN STATEMENT 


'"mw"c£°" 


C 








PRIW 


E NUM 
E (6? 
AT C5 
iHl/l 


BER P 
8 ) ~ 
2H FO 
9X'1.H 


ROBLE 

LLOWI 
2/1 9X 


M 

KI6 IS 
9lH3) 


A LI 


ST OF 


PRIM 


E NUM 


BE 












100 


WRIT 












8 

1 


FORM 
l^'X^ 


•RS 


FROM 


i TO 


1000/ 




101 


I«5 












3 


A-I 


















102 


A = SQ 


RT(A) 






^ J 







u::: 


- 












^il3. 


J=A 


K 


J '2 1 
















104 


DO i 


















105 


1 = 1/ 
















106 


IF(L 


^K-I) 


1?2»4 


























1 


CONT 


INUE 






















JU37 


WRIT 


E (6?1 
At (I 


5)1 






R) 











— 












5 


FORM 


20) 















2 


1 = 1 + 


2 






ERRO 
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IF{1 


000-1 


Y7■)^'> 


L3 












4 


WRIT 


E (6i 


9) 




















9 


FORM 


AT a 


4H PR 


OGRAM 


- - 














-- 




7 


WRIT 


E (6^ 
Af 13 


6) 








i^ THE 








6 


FORM 


IH TH 


is IS 


THE 


END 


PROG 


RAM) 
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STOP 


























END 














-- 
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Figure 2. Sample Program 1 
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SAMPLE PROGRAM 2 

The n points (Xj^, y-^) are to be used to fit an m degree polynomial by 
the least-squares method. 



y = a_+ a^x + aaX^ + ... + an^xi^ 



In order to obtain the coefficients ao/ a^, 
solve the normal equations: 



(1) Woao + W^ai + . 

(2) W^ao + Waai + . 



• + Wpn + najp = Z^ 



, aixi, it is necessary to 



(m+1) Wma + W^ + ia, + 



where: 



• + ^am^m ~ "^m 



W„ = 



w, = 



n 


Z 




f=/^ 


n 






n 


I Xi 

i=1 


Zi 


^ 


i=1 


n 






n 


1 Xi^ 

i=1 


Za 


= 


l^^y,x,= 



Zm = S yiX^m 
i=1 



Warn = ^ Xi^m 
i=1 ' 



After the W's and Z's have been computed, the normal equations are 
solved by the method of elimination which is illustrated by the 
following solution of the normal equations for a second degree polynomi- 
al (m = 2) . 

(1) W^a^ + W^a, + W^aa = Z^, 

(2) W^ao + Waai + Wgaa = Z^ 

(3) Waa^ + W3ai + W^^^a = Za 
The forward solution is as follows: 

1. Divide equation (1) by Wq 

2. Multiply the equation resulting from step 1 by W^ and subtract from 
equation (2) , 

3. Multiply the equation resulting from step 1 by Wa and subtract from 
equation (3) . 
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The resulting equations are: 

(U) a^ + bia^i + ^13^2 = bi* 

(6) 1332^1 + ^3332 = ^34^ 

where: 

baa ~ ^3~^i 3^n r ^33 = W3— busWn , ba^j = Zi^bii^Vif 
baa ~ ''3~t'i a^a # *^33 ~ "i»~bn3Wa # b3«| = Za~b^4jWa 

Steps 1 and 2 are repeated using equations (5) and (6) , with baa ^ind baa 
instead of W© and W^ . The resulting equations are: 

(7) a^ + Casaa = Ca^ 
(o) C33aa ~ C3«t 

where: 

Ca3 ~ has^baa / Ca«» ~ bai»/baa 

C33 = baa-Casbsa # C3«» = bai^-Ca^jbaa 
The backward solution is as follows: 

(9) aa = C3U/C33 from equation (8) 

(10) a^ = Ca«»'"Ca3aa from equation (7) 

(11) a = b^it -bia^i-b^aaa from equation (4) 

Figure 3 is a possible FORTRAN program for carrying out the 
calculations for the case: n = 100, m < 10. Wo» W« # War •••* Warn are 
stored in W (1) , W (2) , W (3) , ..., W(2M+1), respectively. Zq, Z^, Za, 
.-., Zm are stored in Z (1), Z (2), Z (3), ..., Z (M+1) , respectively. 
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FORTRAN Coding Form 
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1 
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Ui^i 
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1 


: 
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DO 5 
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5 


M(J) 
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DO 6 
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■ 
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. 
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1 ■ ■ i ■ 
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Figure 3. Sample Program 2 



IBM 














f ORTRAN Codlno Form 


























.„„r;r;:;. 


'— SAMPLE PROGRAM 2 


rs~fo.5 


G.AfHlC 
















,A0.2 0' 3 


— " l°- G/S6 


PUNC„ 


















"nu™»' J 


FORTRAN STATEMENT 


'"Sl^r" 




7 ! 9 10 11 17 13 H IS 


16' 17 IS 19 20 




73 74 75 76 77 7B 7T N 


16 


W(J) 


= W( 


J)+P 
























17 


DO 2 


1 = 


liLZ 


























DO 2 

J = 


0„X- 
K+I 
I) = 
2 K» 


1?LJ? 
















































20 


BCKi 


W(J-1 


) 
























DO 2 


1»L2 
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Figure 3. Sample Program 2 (Continued) 



Appendix D: Sample Programs 107 



IBM 




































































FOHTRAN Coding Farm 


























































.; 


r>, 


"1 


~ SAMPLE PROGRAM 2 


INSTRLKTIONS 


GRAPHIC 
















PAOE 3 OP 3 


nOCUMMU Idate Qt/Qi^ 


PUNCH 
















CARD ELECTRO NUMIER- 


STATEMENT 
NUMUII 


I 


FORTRAN STATEMENT 


IDENTIFICATION 
SEQUENCE 


"i 2"'3"4'S 




7 8 9 10 11 n 


13 14 


IS 16 17 IS 19 20 21 32 


S3 24 is 26 


i7 28 i9 io 31 32 33 34 


35 36 37 33 39 40 41 42 43 44 49 


11 47 48-49 SO 51 52 S3 54 95 56 57 58 59 60 61 62 63 44 65 66 67 68 69 70 71 72 


73 74 75 76 77 78 79 60 








3 


7 




s 


I 


6 


M 


A 




s 




s 


I 


G 


M 


A 


+ 


B 


( 


I 


- 


1 


•) 


J 


) 


^ 


A 


( 


J 


) 


































































































I 




« 




I 


- 


1 










































































































































A 


( 


I 


) 




s 




6 


( 


I 


t 


L 


B 


) 


- 


S 


I 


6 


M 


A 










































































































W- 







I 


F 




( 


I 


- 


1 


) 




I* 


1 


1 


4 


1 


9 


3 


5 
























- 




- 




















































































4 


1 




^ 


R 


I 


T 


E 




( 


6 


> 


2 


) 




( 


A 


( 


I 


) 


t 


I 


a 


1 


•) 


L 


1 


) 






























































































s 


x 





P 


































- 














































































































E 


N 


D 












































































































































































































































































































































































































































































































































































































































































































































































































































' 




































































































































i 




















































































































































; 





























































































































































































































































































































































































































































































































































































- 


' 






















































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































1 J 3 4 5 


« 


7 e 9 10 11 12 13 14 IS 16 17 le 19 W 21 » » 74 a 24 77 » 29 30 31 33 33 34 IS 34 37 3S 39 « 41 « 43 44 49 4i 47 49 49 so 91 S2 53 S4 55 36 S7 Sa 39 «0 «1 «2 U 64 65 M 67 «e 69 70 71 72 


73 74 75 76 77 78 79 80 1 



^ ttandwd cold fbnn, IBM •laclro B88157, 



Figure 3. Sample Program 2 



set equal to zero. 



L^L and Yj^ are selected. 



The elements of the W array, except W(l), are 
W(l) is set equal to N. For each value of I, X^ 

The powers of Xi are computed and accumulated in the correct W counters. 
The powers of Xi are multiplied by Y^^ and the products are accumulated 
in the correct Z counters. In order to save machine time when the 
object program is being run, the previously computed power of X^ is used 
when computing the next power of Xi. Note the use of variables as index 
parameters. By the time control has passed to statement 17, the 
counters have been set as follows: 



108 



W(l) = N Z(l) = 2 Yj 



N 
2 
1=1 



W(2) = 2 X^. Z(2) = 2 y X^. 
1=1 ^ 1=1 -^ ■ 



N N 

W(3) = 2 Xj^ Z(3) =2 Y^X 2 

1=1 1=1 ^ 



Z(M+1) = 2 Yt-X^-M 
1=1 ^ 



N 
W(2M+1) = 2 X-r^M 
1=1 



By the time control has passed to statement 23, the values of W „ W-l, 
--,, Wam+jL have been placed in the storage locations corresponding to 
columns 1 through M + 1,, rows 1 through M + 1, of the B array, and the 
values of Zo# Zj,» ...» Zm have, been stored in the locations correspond- 
ing to the column of the B array. For example, for the illustrative 
problem (M = 2), columns 1 through 4, rows 1 through 3, of the B array 
would be set to the following computed values: 



Wo 


W:l 


W2 


Zo 


Wi 


Wa 


W3 


Zi 


Wa 


W3 


w,^ 


Z2 


ints 


equat 


ions 


(1). 



This matrix represents equations (l), (2):, and (3), the normal 
equations for M = 2. 

The forward solution, which results in equations (4),, (7), and (8) in 
the illustrative problem,, is carried out by statements 23 through 31. 
By the time control has passed to statement 33, the coefficients of the 
AI terms in the M + 1 equations which would be obtained in hand 
calculations have replaced the contents of the locations corresponding 
to columns 1 through M + 1, rows 1 through M + 1, of the B array, and 
the constants on the right-hand side of the equations have replaced the 
contents of the locations corresponding to column M + 2, rows 1 through 
M + 1, of the B array- For the illustrative problem, columns 1 through 
U, rows 1 through 3, of the B array would be set to the following 
computed values: 



J12 



bi3 bj 







J- »-23 '-'S.lt 

C33 Ca^ 

This matrix represents equations (4), (7), and (8). 

The backward solution, which results in equations (9), (10), and (11) 
in the illustrative problem, is carried out by statements 33 through 40. 
By the time control has passed to statement 41, which prints the values 
of the A9 terms, the values of the (M + 1) *Ai terms have been stored in 
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the M + 1 locations for the A array. For the illustrative problem, the 
A array would contain the following computed values for az, 8i±t and a^, 
respectively: 



Location 


Contents 


A(3) 


C3i»/C33 


A(2) 


Cai»"~C23a2 


A(l) 


bi*f-bi2ai-bi.3a2 



The resulting values of the AI terms are then printed according to 
the FORMAT specification in statement 2,. 
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A format code 59-6 2 
ABS 

(see mathematical function subprograms) 
Addition 

(see arithmetic operators) 
Additional input/output statements 67 

BACKSPACE 67 

END FILE 67 

REWIND 67 
Adjustable dimensions 70-72 
AIMAG 

(see mathematical function subprograms) 
AINT 

(see mathematical function subprograms) 
ALOG 

(see mathematical function subprograms) 
ALOG 10 

(see mathematical function subprograms) 
Alphabetic characters (table) 95 
Alphameric characters 14,95 
AMAXO 

(see mathematical function subprograms) 
AMAXl 

(see mathematical function subprograms) 
American Standards Association (ASA) 

FORTRAN 5 
AMINO 

(see mathematical function subprograms) 
AMINl 

(see mathematical function subprograms) 
AMOD 

(see mathematical function subprograms) 
AND 

(see logical operators) 
Arithmetic and logical assignment 7,27-28 
Arithmetic expressions 19-23 

arithmetic operators 20-23 

mode of 21 

order of computation in 22-23 

use in logical expressions 23 

use of parentheses in 23 
Arithmetic IF 32-33 
Arithmetic operators 20-23 
Arrangement of arrays in storage 19 
Arrays 16-19 

arrangement in storage 19 

declaring the size of 18 

suscripted variables 16 

subscripts 17-18 
ASSIGN statement 31-32 
Assigned GO TO statement 31-32 
ATAN 

(see mathematical function subprograms) 



BACKSPACE statement 67 
Basic input/output statements 

READ Ul-47 

WRITE 47-50 
Basic Operating System 5 
Basic Programming Support 5 
Blank common 76-77 
Blank fields 



40-50 



(see X format code) 
Blank lines 

(see carriage control) 
Blanks 8 
BLOCK DATA subprogram 93 

C (see comments lines) 
CABS 

(see mathematical function subprograms) 
CALL statement 87-88 
Carriage control 66 
CCOS 

(see mathematical function subprograms) 
CDCOS 

(see mathematical function subprograms) 
CDEXP 

(see mathematical function subprograms) 
CDLOG 

(see mathematical function subprograms) 
CDSIN 

(see mathematical function subprograms) 
CDSQRT 

(see mathematical function subprograms) 
CEXP 

(see mathematical function subprograms) 
Character card punch codes 95 
CLOG 

(see mathematical function subprograms) 
CMPLX 

(see mathematical function subprograms) 
Coding form 8 
Comments lines 9-10 
COMMON Statement 74-77 

blank COMMON 76-77 

declaring the size of an array 18 

labeled COMMON 76-77 
Compiler 5 

Complex constants 12-13 
COMPLEX Statement 70-71 

(see also FUNCTION subprograms) 
Computed GO TO statement 30 
CONJG 

(see mathematical function subprograms) 
Constants 7,11-14 

complex 12-13 

double-precision 

integer 10 

literal 13 

logical 13 

real 11-12 
Continuation lines 
CONTINUE statement 
Control statements 

arithmetic IF 



11-12 



7 

37-38 
7,29-39 
32-33 



assigned GO TO 31-32 

computed GO TO 30 

CONTINUE 37-38 

DO 34-36 

END 39 

logical IF 33-34 

PAUSE 38 

STOP 38 
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unconditional GO TO 29 
Conversion codes 

(see format codes) 
COS 

(see mathematical function subprograms) 
CSIN 

(see mathematical function subprograms) 
CSQRT 

(see mathematical function subprograms) 

D decimal exponent 11-12,58-59 

D format code 58-59 

DABS 

(see mathematical function subprograms) 
DATA initialization statement 97-98 
Data set 40 

Data set reference number 40 
DATAN 

(see mathematical function subprograms) 
DATAN 2 

(see mathematical function subprograms) 
DBLE 

(see mathematical function subprograms) 
DCMPLX 

(see mathematical function subprograms) 
DCONJG 

(see mathematical function subprograms) 
DCOS 

(see mathematical function subprograms) 
Decimal exponents 11-12,58-59 
Declaring the size of an array 18 
Device (I/O) 40 
DEXP 

(see mathematical function subprograms) 
DFLOAT 

(see mathematical function subprograms) 
Digit 

(see numeric characters) 
DIM 

(see mathematical function subprograms) 
DIMENSION statement 72-73 

adjustable dimensions 72-73 

declaring the size of an array 18 
Division 

(see arithmetic operators) 
DLOG 

(see mathematical function subprograms) 
DLOGIO 

(see mathematical function subprograms) 
DMAXl 

(see mathematical function subprograms) 
DMINl 

(see mathematical function subprograms) 
DMOD 

(see mathematical function subprograms) 
DO statement 34-37 
DO variable 34-36 
Double-precision constants 11 
DOUBLE PRECISION Statement 98 
DSIGN 

(see mathematical function subprograms) 
DSIN 

(see mathematical function subprograms) 
DSQRT 

(see mathematical function subprograms) 
DTANH 

(see mathematical function subprograms) 
DUMP subprogram 93,102-103 



DVCHK 

(see machine indicator tests) 

E decimal exponent 11-12,58-59 

E format code 58-59 

END FILE statement 67 

END parameter in a READ statement 41 

END statement 39 

in FUNCTION subprograms 85-86 
ENTRY statement 89-91 
EQ (see relational operators) 
EQUIVALENCE Statement 77-79 
ERR parameter in a READ statement 41 
EXIT subprogram 93,102 
EXP 

(see mathematical function subprograms) 
Explicit specification statement 7,70-73 
Exponentiation 23 

(see also arithmetic operators) 
Exponents 

(see decimal exponents) 
Expressions 19-26 

arithmetic 19-23 

logical 23-26 
EXTERNAL Statement 92 

F format code 58 
FALSE 13 

(see also logical expressions) 
Features of operating system FORTRAN IV 

5-6 
Fields 7-8 

blanks (see also X format code) 8 

comments 7-8 

continuation 7 

identification 7 

statement number 7 
FLOAT 

(see mathematical function subprograms) 
Format codes 50-66 

A code 59-61 

carriage control 66 

D and E codes 58-59 

F code 58 

G code 52-57 

H code 62-63 

I code 57-58 

L code 59 

numeric codes 56 

scale factor-P 64-65 

T code 64 

X code 63 
FORMAT statement 40,50-66 

format codes 50-66 

FORTRAN record 40,50-52 

literal data 61-62 

reading FORMAT statements 47 
FORTRAN 

American Standards Association 5 

basic operating system 5 

basic programming support 5 

coding form 7-8 

compiler 5 

library 92,99-103 

object program 5 

record 40,50-52 

source program 5 

statements 7 
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supplied subprograms 92,99-103 
Functions 80-81 

definition of 81 

FUNCTION subprograms 83-86 

reference to 81 

statement function subprograms 81-83 

G format code 52-56 

GE (see relational operators) 

GO TO statements 29-32 

assigned 30-32 

computed 30 

unconditional 29 
GT (see relational operators) 

H format code 62-6 3 

HFIX 

(see mathematical function subprograms) 

Hierarchy of operations 

in a logical expression 25-26 
in an arithmetic expression 22 

I format code 57 
I/O list 

within a NAMELIST 4 3-44,48 

within a READ 41 

within a WRITE 47-4 8 
lABS 

(see mathematical function subprograms) 
IDIM 

(see mathematical function subprograms) 
IDINT 

(see mathematical function subprograms) 
IFIX 

(see mathematical function subprograms) 
Imaginary number 

(see complex constants) 
IMPLICIT specification statement 15,68-70 
In-line 

(see mathematical function subprograms) 
Increment 34 
Indexing I/O lists 46-47 
Indexing parameters in a DO 34-35 
Initial value 34 
Input/output statements 7,40-66 

BACKSPACE 67 

END FILE 67 

READ 41-47 

REWIND 67 

WRITE 47-50 
INT 

(see mathematical function subprograms) 
Integer constants 10 
INTEGER Statements 70-73 

(see also FUNCTION subprograms) 
ISIGN 

(see mathematical function subprograms) 

L format code 59 

Labeled COMMON 7 6-77 

LE (see relational operators) 

Length specification 

(see optional length sp€;cif ication, 
standard length specification) 
Library (see FORTRAN) 
Lines 

(see blank lines, continuation lines) 
List (see I/O list) 



Literal constants 13 

Literal data in a FORMAT statement 61-62 

Logical constants 13 

Logical expressions 23-26 

logical operators 24-25 

order of computation in 25-26 

relational operators 24 

use of parentheses in 26 
Logical IF statement 33-34 
Logical operators 24-25 
LOGICAL statement 70-73 

(see also FUNCTION subprograms) 
Looping (see DO statement) 
LT (see relational operators) 

Machine indicator tests 102 

Mathematical function subprograms 

93,99-101 
MAXO 

(see mathematical function subprograms) 
MAXl 

(see mathematical function subprograms) 
MINO 

(see mathematical function subprograms) 
MINI 

(see mathematical function subprograms) 
Mixed-mode 5-6 

(see also expressions) 
MOD 

(see mathematical function subprograms) 
Mode of an arithmetic expression 19-22 
Muliline listing 56 

Multiple ENTRY into a subprogram 89-91 
Multiplication 

(see arithmetic operators) 

Named common 

(see labeled COMMON) 
NAMELIST statement 42-44,48 
NE (see relational operators) 
Nest of DOS 36-37 
NOT (see logical operators) 
Numeric characters 95 
Numeric FORMAT codes 56-57 

Object program 5 
Operators 

(see arithmetic, logical, relational) 
Optional length specification for 
variables 14-15 

(see also type statements) 
OR (see logical operators) 
Order of computation in arithmetic 

expressions 22-23 
Order of computation in logical 

expressions 25-26 
Out-of-line 

(see mathematical function subprograms) 
OVERFL 

(see machine indicator tests) 

P format code 64 
Parentheses 

use of 23,26 
Parentheses in arithmetic expressions 23 
Parentheses in logical expressions 26 
PAUSE statement '38 
PDUMP subprogram 93,103 
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Pre-defined specification (convention) 15 

PRINT statement 97 

Programming considerations in using a DO 

loop 36-37 
PUNCH statement 96 

Range of a DO statement 3i|-35 
READ statement 41-47 

READ (a) list 45-46 

READ (a,b) list 44-45 

READ (a,x) 42-44 

READ b, list 96 
Reading FORMAT statements 46-47 
REAL 

(see mathematical function subprograms) 
Real constants 11-12 
REAL statement 70-73 

(see also FUNCTION subprograms) 
Relational operators 24 
Repeat constant 43 
RETURN statement 

in a FUNCTION subprogram 84 

in a main program 89 

in a SUBROUTINE subprogram 88-89 
REWIND statement 67 



arithmetic and logical assignment 
7,27-28 

comments lines 8 

continuation of 7 

control 7,29-39 

input / output 7,40-66 

number 7 

specification 7-8,68-79 

subprogram 7,80-93 
STOP statement 38 
SQRT 

(see mathematical function subprograms) 
Subprogram names as arguments 

(see EXTERNAL statement) 
Subprograms 7,80-93 

FORTRAN supplied 92-93,99-103 

FUNCTION 81-83 

statement functions 81-83 

SUBROUTINE 86-89 
Subscripted variable 16-18 
Subscripts 16-18 
Subtraction 

(see arithmetic operators) 
Symbolic unit number 

(see data set reference number) 



Sample program 

program 1 104 

program 2 105-111 
Scale factor-P 64-66 
Sense light subroutines 93 
SIGN 

(see mathematical function subprograms) 
SIN 

(see mathematical function subprograms) 
Slashes in a FORMAT statement 50-51,54-55 
SLITE 

(see machine indicator tests) 
SLITET 

(see machine indicator tests) 
SNGL 

(see mathematical function subprograms) 
Source program 5 
Source program characters 95 
Special characters (table) 95 
Specification statements 7,68-79 

COMMON 74-77 

DIMENSION 72-73 

EQUIVALENCE 77-79 

explicit 70-73 

EXTERNAL 92 

FORMAT 40,50-66 

IMPLICIT 15,68-70 

NAMELIST 42-44,48 
Standard length specification for 
variables 14-15 

(see also type statements) 
Statement numbers 7-8 
Statements 7-9 



T format code 64 
TANH 

(see mathematical function subprograms) 
Test value 35-36 
TRUE 13 

(see also logical expressions) 
Type and length specification 15 
Type declaration 15-16 

explicit 16 

IMPLICIT 15 

pre-defined 15 
Type specification of FUNCTION subprograms 

84-86 
Type statements 68-73 

explicit 70-73 

IMPLICIT 68-70 

Unconditional GO TO statement 29 

Variable FORMAT statements 

(see reading FORMAT statements) 
Variables 7,13-14 

subscripted 16-18 

type declaration 15-16 

types and length specifications 14-15 

variable names 14 

WRITE statement 47-50 
WRITE (a) list 50 
WRITE (a,b) list 49 
WRITE (a,x) 48 

X format code 63 
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